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made it a continuing topic of research for many years. Commercial systems could be 
made more interoperable and military sensors could send demodulated products for 
analysis, to name just two. Noisy channels and multipath fading environments continue 
to make this a challenging problem. This thesis applies classification algorithms that 
have been used in other applications. Nine different digital modulation schemes are 
considered. The criteria for selecting higher-ordered moments and cumulants as features 
for discrimination are discussed. An overview of the classification algorithms considered 
is provided, as well as the statistical models for noisy channels. Results show that the 


scheme proposed here works well in AWGN channels and in moderate fading conditions. 


THIS PAGE INTENTIONALLY LEFT BLANK 


vi 


Il. 


Il. 


IV. 


TABLE OF CONTENTS 


PN TROD UCT ION cis cteciedicesesectctevensidess execast ceed cosetvaialicdesiivesanilunstecieieaieiiess weneiiaetinnets 1 
A. BEC DEV ES costes ssaseerstarecashae cette aa seat sash rican ave aces sata aace ches sk tease stack 1 
B. BACK GROUND ssesiscoie sashes tacts. osscvaceuuiienaéeasasdbeciieusa canst nosasidcoesbiacasdesteeeacnccinatee 1 
C. ORGANIZATION  seccssicotececi sein ccs ces tail acid cccbieti len pee thcetdastecse ce hdeeatibevss ce ttetaee 2 
D. SIMULATION SOP T WARK wiviesscscttecssdotsonscsestdveadsceesdevecdesesesescoseddacdecseisseeses 3 
DIGITAL MODULATION AND WIRELESS CHANNELS ..........ccsccsssssssssseceseees 5 
A. THE NEED FOR MODULATION ..........csccsccsscssccssscscsccccssccsccssccsscsessesssessees 5 
B. DIGITAL MODULATION SCHEMES ...........csscssscsscessseccceccsecssccssscseseneses 6 
1. PVE OIC ENON x oci55eibecoseeshecenssbi cans sucsansoueebenenconsacesscehccmnetovnesacsupeasbesssecenat? 6 
ye Phase) SILC Ie yin 8 oss sccets eciacec ce cawedese dds dacetsswen siecaveevesinndcvessatoveodeceeness 6 
J Quadrature Amplitude Modulation................cccccssssccssssscsssccsescsssessees 8 
4. Frequency-Shift Keyinng............sssccsssssssssssssssssssssssssssssccssesesesssessseseseees 9 
C. SOURCES OF SIGNAL DEGRADATION 2.0... csccssscsccsesceescecssecsscssscsneee 9 
1. FPP OUI CUNO ice csscepscctacaScccareconciss chide seus dnsocaginesdeaguctecesssocseetseesneadiacseaeais 9 
2; Additive White Gaussian Noise ............cssccssssccssssccsssccssscscsessceseessesees 9 
3. Peaelitig Clean S gece ccs caeweccnsives vasuassvacesvsvcncestvesvcnseavectvayaveavsnavevsvcnieens 12 
D. CONCLUSION cuiceisavassenvsctuncsatnsestorgedectes cuvscevstudssouisvonstodeesnssaserasednasubainssnesiao? 14 
SIGNATE: FEA TURES cigscisescisnc ox sisaciccaasnieceiatacaceicbnaden visanibtcassvdeceenveccoisbaatesvieasioracuas 15 
A. IN PIRODU CTO 5cciscoseidecs cates stccanecpicen hiss hac cvsecatcnccicsctlacdheiieieteets epeutantends 15 
B. SPECTRUM BASED CLASSIFICATION OF FSK SCHEMES ............. 15 
C. WOMEN ES ssesésssivcassssitediesnaavs ovis svassseecucessigaessudassucsdcustends siescvecssgecuciasvemdestoserss 16 
D. CUIVIUIGAIN TS ccscscceeciccb ecopce antes viens dcscee eet ctaedctenieeapteeehccconencdi eaves 17 
E. PRA TURE SEEECTION siscccstecsisecesdevesscdencoszetcscacsteteeie aceacusieiavantsdeelacanaens 17 
F. CONCLUSION: ciisatisdseseasccencasancseniastsbscedinetectsscaesstesdesveactsbeaasusesectasepouasinencies 28 
LINEAR AND NONLINEAR CLASSITEFTERG...........cccsccsssssssccssscssccessscsesscsscsssseees 29 
A. INTRODUCTION. sscicsesceecegesicscesh ceccesepes es eenactaeesvses dscisauvccnnveseas sense ececesetsessteuss 29 
B. LINEAR CLASSIFICATION SCHEMES ...........csscsssssssscssscssssssscsesscssseeses 29 
1. Principal Component AmallySis ...........cccssccsssscccssssccssssscsssssssssccsssssees 29 
Zs Linear Discriminant Analysis ............ccsssccssssccsssscssssscssssccssscssssessees 30 
C. NONLINEAR CLASSIFICATION SCHEMES ..........ccscssssssssssssecssccsesseees 32 
1. THe Kernel “Trick” ssscsccevasessves vesszsevxcasvsvosseekvasvenssevnsessvannisssuesvesseesve 32 
Zs Nonlinear Component AnallySis............ccccccsssscsssssccsssccssssscssssssseeceess 33 
3. General Discriminant Amalysis.............csccsssscssssccssssscssssccssssccsssssscees 33 
D. CONCLUSION cisssissessencesstasnctesssenseaissuts sossveshennsvdnsovivlsataadiuasahssenasberssnanstoatas 34 
IMPLEMENTATION AND RESULTS .u.........csccssscsssssssssseccssscssscsssscssescsscsssssssees 35 
A. SIGNAL GENERATION AND CORRUPTION. ............sccscsssssscsscsscssesseees 35 
B. CLASSIFIER IMPLEMENTATION ..........ccssscsssscsssscssscssscsssssesscssscsssssseees 35 
C; LINEAR VERSUS NONLINEAR CLASSIEFTERS...........ccsccssssssssseccsssceees 37 
1. DCA osssaceouicnaccnniasiic iens coansevasinuscanctansediasbanioebundeuas ebanuusndoacevthasvenioaveaaines 38 
2 MUD ss cebeks cass bbivns susene vautohievessbuassoussbaaescusnenunssekeconsstsas souesbansponsettonbabsevese 39 








4. DD Assy cesiesintssoatnwevelscieasiezatd cbiecdal uuidatecauelictousteas he ceavdepmumeamoucbarenmverans 41 
D. PERFORMANCE OF THE PCA AND LDA CLASSIFIERS IN 
AWGN AND FADING CONDITIONS ..........ccccsscssscsssssscscssecsssccscssscssesseees 42 
1. MTV A sos icastiut suahsi sos ccabnbchascascohasesesbeasso cotaat tees shits sousasbostaasennasteatsascbbonsies 42 
2 DOCK vsussesuitisssisidbauss sibcveidodauie iuhiivassassesusieulevsadesivaeeiutduvesdesdidensontcvdaatiere 50 
WE. “CONCTSUSTONS siicspiivucscsusescushhcsivtis seisns ecapituahscsvsunscobanktaaheueseeapi cuskssosevessoauvebsasvavenipne 59 
APPENDIX A: CUMULANTS EXPRESSED AS FUNCTIONS OF EQUAL AND 
LOWER ORDER MOMENT ..........ccccsscsssssssssccssscsssccscssscsscssccsecsscssessssssesssecssscsecees 61 
APPENDIX B: BEHAVIOR OF MOMENTS AND CUMULANTS WITH 
DEC REASTING IONE 5 sisesscnsa issn ctvekscsnucediaiees'ssis iuseiteaby ued sctuaeecssphaviban Sunde cobidueotuinsehessst 63 
A. DWV NONI GY esses snsasaceseaucaa esas cancel nsstetanectatasengsonteceabicassbeestaberoeunttasssmianeasinde 63 
B. AWGN PLUS SLOW, FREQUENCY-FLAT RAYLEIGH FADING......75 
C. AWGN PLUS SLOW, FREQUENCY-FLAT RICEAN FADING ........... 88 
D. AWGN PLUS SLOW, FREQUENCY-SELECTIVE RAYLEIGH 
BEA EDIING wipckicetaecaiadopnads ates sinsctssvedsveessten deoyestiosessstukossesdedieveuducussseviadeustenieavetens 102 
E. AWGN PLUS SLOW, FREQUENCY-SELECTIVE RICEAN 
HA DIING sis pacsnse ss cxntcessaseuseosuienasceasacneshivbeussiedsnsntoastvsoaspsbisedessentesbusoussesssnstesiie 115 
F. AWGN PLUS FAST, FREQUENCY-FLAT RAYLEIGH FADING .....129 
G. AWGN PLUS FAST, FREQUENCY-FLAT RICEAN FADING............ 142 
H. AWGN PLUS’ FAST, FREQUENCY-SELECTIVE RAYLEIGH 
BOA EDIING wo iicaiccieaicacsonsaoiasasinobeavnscdsazesheoncaeuscuckeSuapteaapebavesessv eve eesnass cuss eesareeiaven 156 
I. AWGN PLUS'- FAST, FREQUENCY-SELECTIVE RICEAN 
A DOEING siviasscsdeiswessucassiuss snesienasCeas sunchesnctaavensastensiacessivessocsioausevavovessbdastaievatbonse 169 
APPENDIX C: EFFECT OF CONSTELLATION ROTATION ON MOMENTS 
AND CUMULANTS FOR PSK AND QAM SIGNALS ..........cessssscssscsccsscesceseees 183 
A. IND OUIVEEIIN DS cousSacaccuteh ce dhissusdssbeasietiad se saiceen by vied caseat sdsspaaustassuseticcohicutescutbobiGesocvesdars 183 
1. PS IK SHA isp tcsceseeseerateteseahevencas ease Gre 183 
zs QA WE Siti als siise ss sccsetsssnssseccssvitsensbeacesseontdveseswasenavasucevestaasencepnvesasoveasss 183 
B. IIT EGAN TS sass ciscussecceatiens so sksuseaceaconcesuaessetibcasssobasstsaceasenneomessssuticensbeseseny 184 
1. C,, 4.0 ssssssnsnsnsesnsesnsesncnsesesnsesesssesesnsssnsnsnsesnsesesssesesesessesesesssesesssssssssssnsess 184 
2. C,, 3.1 ssesesesesnsnsesnsnsnsnsnsnsesssnsnsesesnsssesnsnsesnsnsesesesesssesessssssssesssssssssessesess 185 
3. C52 ssesnsnsnsnsnsnsnsnsnsnsnsnsnssnsesesnsnssesesessesesssesesesesessesssesesssssssssssesnsnss 185 
4. Cy, gp aasavsnraccsernrdacesssnsacesdssncnosdurntenssivensssersscdoeentuadesennscseswesdsqessvonseusebens 185 
5. CC, 5,1 coessecsssssscncecnssccecssssscssssscosesnsnsceessesscsssssnsosssnsssccssesecsssssnsesesncsecnees 185 
6. Cy, 4.2 ssssnsnsnsnsnsnsnsnsnsnsssnssnsesesesessnsesnsesnsesnsesesesesesesesssssssesesssssssesssnsess 185 
7. Cy, 4.3 sssesnsesesnsnsnsnsnsnssnsnsnsesnsnsesesnsesnsesnsnsnsesasesssessessesssssesesesssssssssssess 186 
8. C,, 9.0 ssnsesnsesesesesnsnsnsesesesesesesesesesesesesesesesesesssesesesesesssesesssssesssssssssssesnsnss 186 
9. Cy 71> Cy6.2» C531 ANA Cy 4.4 sorscscsssccecsrscsccscnsesnssesssscesersssscscnsesessassssees 186 
APPENDEX D? MATLAB CODES isscceciccsceeatuacsusscvesdieatwcseuscesesciocaedinessouoessouseciieesevescenseeses 187 


IDINEAR CLASSIPIERS wéssescoctsctas sssvieasecuvedeveecevedeasstansidctdsseuscepstsdscasoucadetis 187 


1. DU @SCr ipl Onis 5scadssevesccsceveccheceevccs Geccacpy sees bcs ieeaivtesacnnahccaventecedeavidass 187 
2 BC Ara sca cneehss satan tasiac ae cscudabentanskala ssaeaebtesebeabeanaacaeissacasantasasiTeaemensione 187 
3 DADA SW ss ayissieshcasiwanbonisiisevandwibeaswnashousceassovasstiesesavicasendteasebebonsstuassedepuese 188 
4 LINEAR TES FER sUitscsccéoscicsasvensssesspestocatoncabosiboreibseppsovsanvansvassnnseseaty 190 
NONLINEAR CLASSIFIER ........c:csccsssssossssccsscsccssccsesssossescsessescoessecsssesses 194 
1. TV@SCETPEIGTIS  sace5 coisas ceecaxevcccucescstesavouseveecdenecuwebeaveuneorsunaceudesieeeuaceuennens 194 
2. FBC Asi gsssciecsasvcusvcsptiveksctnseueas svasesantvartecabhcunkosvesessasnyeovtdesocns tasprensians 195 
3. UDA INT cus isseh i cisiceasagsievannd cnc tbesasiaeiaubstaassatiasenstodocabessaucadpecduadsesbaaeeasponce 196 
4. NONEEN TES DE Rin scsicasictssceasiencabisesescastansexcaseabceabivesvcusiibesagsscons 198 
SUPPORTING FUNC TIONS ws ssesisevsttevssdondsonstivscvbadestasdesvstesssdescsounstsosdeveadone 203 
L. DGSCEEDUION sds ssccvicscesicsces dusnacdcusitvacdussissxdesesesnssdeieavsecssnadssdenerns satvaiaadiees 203 
2 PRO EG ENGI gsicascsssncsiscanh cual scucsesnssetcsseacsouesecspocwes sehusobecrobibebanseuioecte 204 
3 PICAENPROE'S sii sy sasvassatcscensconwcccdessucadpeatocdoestveceesssacadpeadenssenaaseannoncs 206 
4 FREER NED ai iiss seca deatssihdsvaticchessuniesbeaniasieacisiiteustheassnsteadebeneaaceasaouooee 208 
5 DV TDD ES En si 6sscs sh caecedessekcepsccesasteassecnss cnstelesseubscussusssbrussassess conven saute thinks 208 
LIST OF REFERENCES sieisssassscvedessverssesasonessseivecsssousstesssenscasusevassobecvecsdesssetespouvscosnvenasensas 211 
ENEDIAL DES ERB UTIONIIGINS D sicosiicesscccecascasecastnesdaes cathecbendeeedeneetocaashebeatandetdevsensetaenteeaks 213 


1X 


THIS PAGE INTENTIONALLY LEFT BLANK 


Figure 1. 
Figure 2. 
Figure 3. 
Figure 4. 
Figure 5. 
Figure 6. 
Figure 7. 


Figure 8. 
Figure 9. 


Figure 10. 
Figure 11. 
Figure 12. 
Figure 13. 
Figure 14. 
Figure 15. 


Figure 16. 
Figure 17. 


Figure 18. 
Figure 19. 
Figure 20. 
Figure 21. 
Figure 22. 
Figure 23. 
Figure 24. 
Figure 25. 
Figure 26. 
Figure 27. 
Figure 28. 
Figure 29. 
Figure 30. 
Figure 31. 
Figure 32. 


LIST OF FIGURES 


OPSK  Coriste Nations 5 ssssccssercaivsisssiasieseechishesed sasseeeealaasisneorerenla oct aa teieeees 7 
16-QAM Constellation. ...........cccccccssssssccececececsesesseaeceseceesesesssaeeececceenessnseaeeesees 8 
Power Spectral Density Of AWGN. uu... eeeeeecssceceseceesscceenscceensececnsceeeecens 10 
Effect of AWGN on 16-QAM (SNR = 20dB). 0... eeeesesseeereeeeceteceaeeneeeneeneeens 11 
Effect of Rayleigh Fading on 16-QAM (SNR = 20dB)....... eee eeeeeeeeeseeeeneeees 13 
Single-Sided Spectrum of BFSK. .0..... cece ceececsececseeeceeeeeceeeeeceeeeeeseeeenaeeeenas 16 
Cg KAWGIN Only) sas aA ecole Ae ee we ee eects 23 
Ce CS WGN, OM) na vant se eed ze et i vasdaa alk teal catcada eral esemet tania eee 24 
Gg WGN Oy) nto cd teen res hea aa eie eae aera acim tastes 24 
C 4g (slow; flat Rayleigh tagging) s.icestuce.t acesesttartec seat riaatiaherk teint 25 
C24 (slow, flat Rayleigh fading ) x ccscac. wigs anti eeescpsuncacuausecscniatesetnshstatadeibennn 20 
C.94 Glow, flat Rayleigh fading )si3:.:cccehse.cresedeaethacssiemerienncienjeeceieeesiniee 26 
C.49 (fast, frequency-selective Rayleigh fading). ........:cceseeeeeeeseeeeeees 26 
C,;, (fast, frequency-selective Rayleigh fading).........cscsseeeeeeeeeeeeeees ZL 
C,.¢9 (fast, frequency-selective Rayleigh fading). .........:scsseeseeeeeseeeeeees 2] 
Classification Flow Chart for Test Signals .0....... eee eeeeesseceesteeeesteeeeneeeeaeees a1 
Ee STW GIN aS sete Ra aa AE a Rea 63 
EE ie MIL WIGIN, hse etesd xasscoehvchcaraannest peaiesaak nie acl ee eee loo 63 
IE ys A A WAIN ot scabs cen push euehcacaagavoel eae eam aioe dickens 64 
IEE ch, WSR NV CIN ce cse nei ged cstsctasn sy edtace nad ea ati ead aauen ORT ona eae potted Ue 64 
Ey ANT PN WGN acct 222s aie taader a Pelee caete vest tig Saea psa ret eed a eee cated rae 65 
Fogg: MDW GIN, 52 ccsist ects cosdotertaeceadiap erriedeeeeianaeedeueaiean) ead he caseeeaee: 65 
Ege TUN WIG, ace ecsigeccgeneea tbadiisnaaiutese eGR ca adiean amuse ceviie te teak tetas uaneaetmtaneones 66 
Bi go: AW OW RING Scacatiss dices SSaseecnsinnes ace joeanty sevesessaacyaeenesdaciasa Sueatepsadeensaereeeiatans 66 
t Diperag eve Gt Enero cor teeeere reece ee tre etree ace see reer earner 67 
ME og NTLZEN WW GOIN ech tect store FI cng vedas Agana epc RA a SA Mn SS 67 
MULE WN GIN hci ee Pat Aaa PON cate ata 68 
Bo getgs AIM PW Ct Nos ti hea inectlbate tunhadanttor vant ew adres gine pars tunghes dawn ouvacziaaned Simetipare tecahestae 68 
EE 8 ATV PR WWIGIIN ioe as cain eas ga seats dana GA casas aes aU MRe Naa eee ea 69 
see MUN WGN tite i catch t esta canta ee acca cet eel eta ee ko etaeagetths 69 
Ce pfoM WGN tata tiated hat cae eo atric Sai ashi hap Jae ateaced nen ctorer, 70 
Cy AT PA WIG os exh tres tavingrag edn parceneaa tay eass Ra taa opine a eatinan ina avog Riana Goda deeadye 70 


Xl 


Figure 33. 
Figure 34. 
Figure 35. 
Figure 36. 
Figure 37. 
Figure 38. 
Figure 39. 
Figure 40. 
Figure 41. 
Figure 42. 
Figure 43. 
Figure 44. 
Figure 45. 
Figure 46. 
Figure 47. 
Figure 48. 
Figure 49. 
Figure 50. 
Figure 51. 
Figure 52. 
Figure 53. 
Figure 54. 
Figure 55. 
Figure 56. 
Figure 57. 
Figure 58. 
Figure 59. 
Figure 60. 
Figure 61. 
Figure 62. 
Figure 63. 
Figure 64. 
Figure 65. 


Cog IN AWGIN, 2 aye eee Sakae cae eee as Sela Geese eas 71 
Gag MUL WGN pce et rasan datos toanhoanatatreesoa uta nasa teneeaa ata 71 
C25 TN DW IGIN Go sesaset hc adi dp oat stesso siete aries oad ed SL aces pa cdatees de 
Co ggy MD AWGN cc gees sietteateeatued tiatanisautieni ine toa tapeeltendie chzeivaneubtaneataananeendeidmotate 12 
Gam TOWN, .sccccttates chatesaratiisi attend sceeental acta nseeaaleectatereectinteciasseaiuiaice 73 
Oop AWG Ne sect ciatac andes tata iat ete alaced a tarelibs acca as tataiee adsl eases 73 
es AWGN cde os tarde adel ae, Se tea a at ceded sa eas a ea So iadaaad anda Doel daeaa 74 
Copy. ATA WGN, sesssveheiaanospselgrahertiayhasmige iments tenet eatin Hen ORAS OIE 74 
gp UM PR WGN acct senporah nel aes tlensk natant cstagdulned daahadeaketanlaracnraanninelaralaoeu 15 
dg: MN AWGIN casi Sas tae naa ahied eee aaa G ars Sens ele = 
E,,) in AWGN and Slow, Frequency-Flat Rayleigh Fading..............0 76 
E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............:: 76 
E49 11 AWGN and Slow, Frequency-Flat Rayleigh Fading.............:: 58 
E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. .............. i) 
E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading.............:8 78 
E69 11 AWGN and Slow, Frequency-Flat Rayleigh Fading..............0 78 
E.;, in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............: 79 
E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading.............. 79 
E,,; in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............:: 80 
E...9 in AWGN and Slow, Frequency-Flat Rayleigh Fading...............: 80 
E,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............:: 81 
E.¢. in AWGN and Slow, Frequency-Flat Rayleigh Fading.............. 81 
E.;; 1 AWGN and Slow, Frequency-Flat Rayleigh Fading.............: 82 
E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading.............:.0 82 
C.49 in AWGN and Slow, Frequency-Flat Rayleigh Fading.............0 83 
C,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............. 83 
C,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading............:: 84 
C.69 11 AWGN and Slow, Frequency-Flat Rayleigh Fading.............. 84 
C,;, in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............:: 85 
C,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading.............: 85 
C.,; in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............. 86 
C..9 in AWGN and Slow, Frequency-Flat Rayleigh Fading... 86 
C,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. ............:c 87 


Xii 


Figure 66. 
Figure 67. 
Figure 68. 
Figure 69. 
Figure 70. 
Figure 71. 
Figure 72. 
Figure 73. 
Figure 74. 
Figure 75. 
Figure 76. 
Figure 77. 
Figure 78. 
Figure 79. 
Figure 80. 
Figure 81. 
Figure 82. 
Figure 83. 
Figure 84. 
Figure 85. 
Figure 86. 
Figure 87. 
Figure 88. 
Figure 89. 
Figure 90. 
Figure 91. 
Figure 92. 
Figure 93. 
Figure 94. 
Figure 95. 
Figure 96. 
Figure 97. 
Figure 98. 


mom mA Af 


a 


I oe oe a 


Ya A 


Pe bah Ol CY ea Ey, Oy Oey sO Cy 


Sie 


62 In AWGN and Slow, Frequency-Flat Rayleigh Fading.............. 87 
153 12 AWGN and Slow, Frequency-Flat Rayleigh Fading..............:+ 88 
44 In AWGN and Slow, Frequency-Flat Rayleigh Fading..............:00 88 
in AWGN and Slow, Frequency-Flat Ricean Fading. «0.0.0... eeeeeeeee 89 

i, 1n AWGN and Slow, Frequency-Flat Ricean Fading..........:.:ceee 90 
in AWGN and Slow, Frequency-Flat Ricean Fading. ......... cesses 90 

3, In AWGN and Slow, Frequency-Flat Ricean Fading. ..........: see 91 
3. In AWGN and Slow, Frequency-Flat Ricean Fading. ..........: scree 91 
in AWGN and Slow, Frequency-Flat Ricean Fading. 0.0.0... eeeeeeeee 92 

, in AWGN and Slow, Frequency-Flat Ricean Fading. ...........:ceeeee 92 
in AWGN and Slow, Frequency-Flat Ricean Fading. «0.0.0... eeeeeeeee 93 

in AWGN and Slow, Frequency-Flat Ricean Fading. ...........e ee eeeeeeee 93 

in AWGN and Slow, Frequency-Flat Ricean Fading. ......... cesses 94 

, in AWGN and Slow, Frequency-Flat Ricean Fading. ..........: cee 94 
in AWGN and Slow, Frequency-Flat Ricean Fading. ......... eee 95 

in AWGN and Slow, Frequency-Flat Ricean Fading. 0.0.0... eee 95 

in AWGN and Slow, Frequency-Flat Ricean Fading. 0.0... cece 96 

in AWGN and Slow, Frequency-Flat Ricean Fading. 0.0... eeeeeeeee 96 

3, 1 AWGN and Slow, Frequency-Flat Ricean Fading. .........:.:cceeeceee 9] 
x22 19 AWGN and Slow, Frequency-Flat Ricean Fading. ......... sce aT 
690 1 AWGN and Slow, Frequency-Flat Ricean Fading. «0.0... 98 
51 11 AWGN and Slow, Frequency-Flat Ricean Fading. ...........ceeeees 98 
42 in AWGN and Slow, Frequency-Flat Ricean Fading. «0.0... 99 
133 In AWGN and Slow, Frequency-Flat Ricean Fading. .........ceeees 99 
in AWGN and Slow, Frequency-Flat Ricean Fading. ......... eee 100 

‘7, 19 AWGN and Slow, Frequency-Flat Ricean Fading. .......... eee 100 
62 11 AWGN and Slow, Frequency-Flat Ricean Fading. ..........::eeee 101 
53 In AWGN and Slow, Frequency-Flat Ricean Fading... 101 
44 In AWGN and Slow, Frequency-Flat Ricean Fading. ...........::eeee 102 
29 in AWGN and Slow, Frequency-Selective Rayleigh Fading.............. 103 
t11 12 AWGN and Slow, Frequency-Selective Rayleigh Fading............... 103 
4o In AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 104 
3, In AWGN and Slow, Frequency-Selective Rayleigh Fading............... 104 


Xlll 


Figure 99. 


Figure 100. 
Figure 101. 
Figure 102. 
Figure 103. 
Figure 104. 
Figure 105. 
Figure 106. 
Figure 107. 
Figure 108. 
Figure 109. 
Figure 110. 
Figure 111. 
Figure 112. 
Figure 113. 
Figure 114. 
Figure 115. 
Figure 116. 
Figure 117. 
Figure 118. 
Figure 119. 
Figure 120. 
Figure 121. 
Figure 122. 
Figure 123. 
Figure 124. 
Figure 125. 
Figure 126. 
Figure 127. 
Figure 128. 
Figure 129. 
Figure 130. 
Figure 131. 


ae 
N 


oY 
° 


= om om 


= 
N 


by ty 


& 
we 
es) 


ty 


a N 
ho 


nn 
Ww 


care ese SS 


ea 
aN 


eS 
a 


a 
& 
~ 


Ae) 
b 
8 


a 


Sa 
° 


aS 
w 
2 


Ae 
nS 
& 


io) 


* 


QD 


= 


ao 


oy 
re) 


ee 
Sh 
ne 


12) 
" 
c 


DS 
° 


by om 


by 


=m 


ty 
Bb S) 
N in N 


we 
w 


Bae a 


ty 


in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 105 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 105 
in AWGN and Slow, Frequency-Selective Rayleigh Fading............... 106 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 106 
in AWGN and Slow, Frequency-Selective Rayleigh Fading............... 107 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 107 
in AWGN and Slow, Frequency-Selective Rayleigh Fading............... 108 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 108 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 109 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 109 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 110 
in AWGN and Slow, Frequency-Selective Rayleigh Fading............... 110 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 111 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 111 
in AWGN and Slow, Frequency-Selective Rayleigh Fading............... 112 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 112 
3 in AWGN and Slow, Frequency-Selective Rayleigh Fading............... 113 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 113 
, 1n AWGN and Slow, Frequency-Selective Rayleigh Fading............... 114 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 114 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 115 
in AWGN and Slow, Frequency-Selective Rayleigh Fading. ............. 115 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 116 
1, 1n AWGN and Slow, Frequency-Selective Ricean Fading. ................ 117 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 117 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 118 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 118 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 119 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 119 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 120 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 120 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 12 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 121 


XIV 


Figure 132. 
Figure 133. 
Figure 134. 
Figure 135. 
Figure 136. 
Figure 137. 
Figure 138. 
Figure 139. 
Figure 140. 
Figure 141. 
Figure 142. 
Figure 143. 
Figure 144, 
Figure 145. 
Figure 146. 
Figure 147, 
Figure 148. 
Figure 149, 
Figure 150. 
Figure 151. 
Figure 152. 
Figure 153. 
Figure 154. 
Figure 155. 
Figure 156. 
Figure 157. 
Figure 158. 
Figure 159. 
Figure 160. 
Figure 161. 
Figure 162. 
Figure 163. 
Figure 164. 


=m om 


Ya 


= 


Q 


ca 


Q 


ca 


io) 


jo) ne ca 


= 


io) 


a 


RGo Slew Dele an 


ty 


colossal 


by 


om om om 


io) 


Cds 


a 


62 In AWGN and Slow, Frequency-Selective Ricean Fading.................. 122 
53 In AWGN and Slow, Frequency-Selective Ricean Fading.................. 122 
44 in AWGN and Slow, Frequency-Selective Ricean Fading.................. 123 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 123 

3, In AWGN and Slow, Frequency-Selective Ricean Fading.................. 124 
32 In AWGN and Slow, Frequency-Selective Ricean Fading.................. 124 
«9 1n AWGN and Slow, Frequency-Selective Ricean Fading.................. 125 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 125 

4. in AWGN and Slow, Frequency-Selective Ricean Fading.................. 126 
3 in AWGN and Slow, Frequency-Selective Ricean Fading.................. 126 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 127 

, 1n AWGN and Slow, Frequency-Selective Ricean Fading.................. 127 
in AWGN and Slow, Frequency-Selective Ricean Fading.................. 128 

in AWGN and Slow, Frequency-Selective Ricean Fading.................. 128 

in AWGN and Slow, Frequency-Selective Ricean Fading.................. 129 

in AWGN and Fast, Frequency-Flat Rayleigh Fading. ............. ee 130 

1, 1n AWGN and Fast, Frequency-Flat Rayleigh Fading... 130 
in AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. 131 

3, In AWGN and Fast, Frequency-Flat Rayleigh Fading.............cee 131 
3. In AWGN and Fast, Frequency-Flat Rayleigh Fading. .............:+ 132 
so 1n AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. 132 
5, in AWGN and Fast, Frequency-Flat Rayleigh Fading.............e 133 
42 in AWGN and Fast, Frequency-Flat Rayleigh Fading. .............06 E33 
33 In AWGN and Fast, Frequency-Flat Rayleigh Fading.............cce 134 
in AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. ee 134 

7, In AWGN and Fast, Frequency-Flat Rayleigh Fading.............. 135 
62 In AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. 135 
53 In AWGN and Fast, Frequency-Flat Rayleigh Fading. .............06 136 
44 in AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. 136 
in AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. ee 137 

, In AWGN and Fast, Frequency-Flat Rayleigh Fading.............. 137 
‘x22 19 AWGN and Fast, Frequency-Flat Rayleigh Fading. ..............:+ 138 
so 1n AWGN and Fast, Frequency-Flat Rayleigh Fading. ...........:. 138 


XV 


Figure 165. 
Figure 166. 
Figure 167. 
Figure 168. 
Figure 169. 
Figure 170. 
Figure 171. 
Figure 172. 
Figure 173. 
Figure 174. 
Figure 175. 
Figure 176. 
Figure 177. 
Figure 178. 
Figure 179. 
Figure 180. 
Figure 181. 
Figure 182. 
Figure 183. 
Figure 184. 
Figure 185. 
Figure 186. 
Figure 187. 
Figure 188. 
Figure 189. 
Figure 190. 
Figure 191. 
Figure 192. 
Figure 193. 
Figure 194, 
Figure 195. 
Figure 196. 
Figure 197, 


one) 


= 


Q 


oem 


ey 
Na 
a 


Q 


= 


Q 


ca 


io) 


Ba 


ty 


gi Sy EN eh er 


Oo mm mm 


io) 


= 


oe) 


cad 


Q 


= 


Q 


Q ne ad 


oun om 


ca 


Q 


8 
a 
re) 


A 


ca 


in AWGN and Fast, Frequency-Flat Rayleigh Fading.........00.. eee 139 
42 In AWGN and Fast, Frequency-Flat Rayleigh Fading. ............. 139 
33 In AWGN and Fast, Frequency-Flat Rayleigh Fading..............::0 140 

in AWGN and Fast, Frequency-Flat Rayleigh Fading. ............. ee 140 

in AWGN and Fast, Frequency-Flat Rayleigh Fading.........0.. eee 141 
62 In AWGN and Fast, Frequency-Flat Rayleigh Fading. ...............0 141 
53 in AWGN and Fast, Frequency-Flat Rayleigh Fading. .............. 142 
44 In AWGN and Fast, Frequency-Flat Rayleigh Fading. ............ 142 

in AWGN and Fast, Frequency-Flat Ricean Fading...........eeeeeeeeees 143 
i, 1n AWGN and Fast, Frequency-Flat Ricean Fading. ..........::ceeeee 144 
49 in AWGN and Fast, Frequency-Flat Ricean Fading.......... cee 144 
3, In AWGN and Fast, Frequency-Flat Ricean Fading. ..........:cceeeeee 145 
32 In AWGN and Fast, Frequency-Flat Ricean Fading...........:cceeee 145 
59 1n AWGN and Fast, Frequency-Flat Ricean Fading............eeeeee 146 
5, in AWGN and Fast, Frequency-Flat Ricean Fading... 146 
4. 1 AWGN and Fast, Frequency-Flat Ricean Fading...........:cceeee 147 
33 In AWGN and Fast, Frequency-Flat Ricean Fading. .........:cceeeeeee 147 

in AWGN and Fast, Frequency-Flat Ricean Fading...........ceeeeeeeeees 148 
7, in AWGN and Fast, Frequency-Flat Ricean Fading. .........:cceeee 148 
6. In AWGN and Fast, Frequency-Flat Ricean Fading.............sseee 149 
53 In AWGN and Fast, Frequency-Flat Ricean Fading...........cceeees 149 
44 in AWGN and Fast, Frequency-Flat Ricean Fading... 150 
49 in AWGN and Fast, Frequency-Flat Ricean Fading... 150 
3, In AWGN and Fast, Frequency-Flat Ricean Fading. ..........::ceecee 151 
22 in AWGN and Fast, Frequency-Flat Ricean Fading............cssceeee 151 
60 1n AWGN and Fast, Frequency-Flat Ricean Fading........... eee 152 
5, in AWGN and Fast, Frequency-Flat Ricean Fading... 152 
4. in AWGN and Fast, Frequency-Flat Ricean Fading... 153 
33 In AWGN and Fast, Frequency-Flat Ricean Fading. ..........:ccceseeee 153 

in AWGN and Fast, Frequency-Flat Ricean Fading...........eeeeeeeees 154 
;, In AWGN and Fast, Frequency-Flat Ricean Fading...........:eeeeee 154 

in AWGN and Fast, Frequency-Flat Ricean Fading........... eee eeeeeeees 155 
53 In AWGN and Fast, Frequency-Flat Ricean Fading... 155 


XVi 


Figure 198. 
Figure 199, 
Figure 200. 
Figure 201. 
Figure 202. 
Figure 203. 
Figure 204. 
Figure 205. 
Figure 206. 
Figure 207. 
Figure 208. 
Figure 209. 
Figure 210. 
Figure 211. 
Figure 212. 
Figure 213. 
Figure 214. 
Figure 215. 
Figure 216. 
Figure 217. 
Figure 218. 
Figure 219. 
Figure 220. 
Figure 221. 
Figure 222. 
Figure 223. 
Figure 224. 
Figure 225. 
Figure 226. 
Figure 227. 
Figure 228. 
Figure 229, 
Figure 230. 


2 


a 
-s 
mn 


D. 
° 


by om 


ty 


x,4,0 


BF BR Dp RK Bw 
No = oO NH 


om om om om om 


we 
ioe) 


by 


x,8,0 


Ny 
2S. 


3543 


et 


Q 


fad 


AoA 


3252 


af 


io) 


a 


Sl 


oe) 


ca 


04,2 


Q 


& 


Hes) 


Qa 


= 


io) 


ca 


6,2 


Q 


a 


io) 


a 


by 


x,2,0 


Sines 


POEL. 


Se: a 


oy 
° 


in AWGN and Fast, Frequency-Flat Ricean Fading...........ceeeeeeeeeees 156 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 157 
1, 1n AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 157 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 158 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 158 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 159 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 159 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 160 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 160 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 161 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 161 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 162 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 162 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 163 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 163 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 164 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 164 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 165 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 165 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 166 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 166 
in AWGN and Fast, Frequency-Selective Rayleigh Fading. ............... 167 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 167 
, 1n AWGN and Fast, Frequency-Selective Rayleigh Fading................ 168 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 168 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 169 
in AWGN and Fast, Frequency-Selective Rayleigh Fading................ 169 
in AWGN and Fast, Frequency-Selective Ricean Fading. ............... 170 
1, 1n AWGN and Fast, Frequency-Selective Ricean Fading...............0++ 170 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................ 171 
in AWGN and Fast, Frequency-Selective Ricean Fading. .................. 171 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................. 172 
in AWGN and Fast, Frequency-Selective Ricean Fading. .................. 172 


XVii 


Figure 231. 
Figure 232. 
Figure 233. 
Figure 234. 
Figure 235. 
Figure 236. 
Figure 237. 
Figure 238. 
Figure 239, 
Figure 240. 
Figure 241. 
Figure 242. 
Figure 243, 
Figure 244, 
Figure 245. 
Figure 246. 
Figure 247. 
Figure 248. 
Figure 249, 
Figure 250. 


by 


by om 


we 
te 


nn 
Ww 


mI Hoe a 


ae 
aN 


Qa 


tad 
i) 
— 


ee 
b 
b 


Ya 


bad 
wn 
— 


a 


- 
ty 


OA 
bw 
eo 


one) 


# 
Na 
7 


ba 
Sh 
N 


2 


2 


& 
- 
aN 


in AWGN and Fast, Frequency-Selective Ricean Fading. .............. 173 
in AWGN and Fast, Frequency-Selective Ricean Fading. .................. 173 
in AWGN and Fast, Frequency-Selective Ricean Fading. .............. 174 
in AWGN and Fast, Frequency-Selective Ricean Fading. ............0.. 174 
, In AWGN and Fast, Frequency-Selective Ricean Fading. ..............+ 175 
in AWGN and Fast, Frequency-Selective Ricean Fading. ............0... 175 
in AWGN and Fast, Frequency-Selective Ricean Fading. ............... 176 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................ 176 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................ 177 
in AWGN and Fast, Frequency-Selective Ricean Fading. .............. 177 
in AWGN and Fast, Frequency-Selective Ricean Fading. .................. 178 
in AWGN and Fast, Frequency-Selective Ricean Fading. .............. 178 
in AWGN and Fast, Frequency-Selective Ricean Fading. .............. 179 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................ 179 
in AWGN and Fast, Frequency-Selective Ricean Fading. .............. 180 
in AWGN and Fast, Frequency-Selective Ricean Fading. .................. 180 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................ 181 
in AWGN and Fast, Frequency-Selective Ricean Fading. ................ 181 
in AWGN and Fast, Frequency-Selective Ricean Fading. .............0 182 
in AWGN and Fast, Frequency-Selective Ricean Fading. ............... 182 


XVill 


Table 1. 
Table 2. 
Table 3. 
Table 4. 
Table 5. 
Table 6. 
Table 7. 
Table 8. 
Table 9. 


Table 10. 


Table 11. 


Table 12. 


Table 13. 


Table 14. 


Table 15. 


Table 16. 


Table 17. 


Table 18. 


Table 19. 


Table 20. 


Table 21. 


Table 22. 


Table 23. 


Table 24. 


Table 25. 


Table 26. 


LIST OF TABLES 


Moments of Digitally Modulated Signals (SNR=©0 ). .... eee eeseeeeeseeenteeeteeeees 19 
Cumulants of Digitally Modulated Signals (SNR=0© )..... ee eeeeseeeseeeeeeeees 20 
Normalized Cumulants of Digitally Modulated Signals (SNR=00).............. 21 
Confusion Matrix for PCA classifier in AWGN, SNR = 20 dB... 38 
Confusion Matrix for PCA classifier in AWGN, SNR =5 dB... 39 
Confusion Matrix for LDA classifier in AWGN, SNR = 20 dB. uo... 39 
Confusion Matrix for LDA classifier in AWGN, SNR =5 dB. woe 40 
Confusion matrix for KPCA classifier in AWGN, SNR = 20 dB. ........ 40 
Confusion matrix for KPCA classifier in AWGN, SNR =5 dB. uu... 41 
Confusion Matrix for GDA classifier in AWGN, SNR = 20 dB... 41 
Confusion Matrix for GDA classifier in AWGN, SNR =5 dB... 42 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Flat Rayleigh Fading: (SNR = 20 dB), sciiccescesacesascedesdancasacncacecaccelssincwsscadarsenaye 43 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Blat:Rayleigh Padding (SNR =Ss0B), 2 iscic hd Ghee ia teenie ke 43 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Flat Ricean Fading (SNR = 20 dB)... eeeeeesccecsseceeeseeceseeeceeeeeeseeesseeeenaeeees 44 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Plat Ricean Fading (SNR: = 3 dB) icc2ieccseas deseenieeectendd A eieea tenes: 44 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleigh Fading (SNK = 2 OCdB ) eset octncbecias eects vvgueleaeupakcncumestlavatehs beacuconen’ 45 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleich Padine (SINR = 5, dB), 23 ach access sceaaaicauctisdicacaneaeteotes aacangeeeennvase 45 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat 
Ricean: Fading (SNR = 20 GB yy s2cnnicbor ele let oeal Resta ati paeceuictoia uae 46 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat 
Riceart Paint SINR 5 CIB 3 ics ped ccc sssaepcna let dacdusilauetei egnaceedelessiivoe be stnosse: 46 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Selective Rayleigh Fading (SNR = 20 CB)... eee eeeeeeeseecsneceeeeeneeenaeenseeneees 47 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Selective Rayleigh Fading (SNR = 5 dB). ......e eee eeeceeeeseceseeeecsneeeeeneeeeneeeenes 47 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Selective Ricean Fading (SNR = 20 dB). ..0......e eee eeeceeccesoeeesseeteeceeaeesoneesseeeees 48 
Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency- 
Selective Ricean Fading (SNR = 5 dB). 0... eeeceeescecssececsecceeeseeeeseeeeeseceenaeeesaes 48 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency- 
Selective Rayleigh Fading (SNR = 20 CB)... eee eeeeeeeseeceneceeeeeneecnneenseensees 49 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency- 
Selective Rayleigh Fading (SNR = 5 dB)... cece eeeeeeceeeseecneeeseeeeseecnaeeneeeenes 49 
Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency- 
Selective Ricean Fading (SNR = 20 GB). .......cccccecseccessceceeeeeenseecesteeseneesenteeee 50 


XIX 


Table 27. 


Table 28. 


Table 29. 


Table 30. 


Table 31. 


Table 32. 


Table 33. 


Table 34. 


Table 35. 


Table 36. 


Table 37. 


Table 38. 


Table 39. 


Table 40. 


Table 41. 


Table 42. 


Table 43. 


Table 44. 


Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency- 
Selective Ricean Fading (SNR = 5. dB), .........sccsssccsssescessnsescsescessceseneeceenstonsnes 50 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Plat Rayleigh Fading (SNR = 20 dB). sic. 4.ccn cdi eee eit 51 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Flat Rayleigh Fading (SNR = 5 CB). 0... eee eeeeececsseceeeeeeceeeeeceeeeeceeeesseeeenaeeees 51 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Plat: Ricean Fading (SNR =20'0B ) 3 tisssot xsgeonitidaniegdipicetield asenns a2 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Blat-Ricean Padinge (SNR = 5 GB) siiccsiavceessdesdecedeacdstenscavavsonsuances eatencestiavatoonbes 52 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleigh Fading (SNR:=20 CB), 4. atssccce) sasdoeaaics gusees sabsndandecerceusaqaceteassetebetetene ae, 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleigh Fading (SNK = 5 dB) jew ice een eid alee aie eae 53 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat 
RiceanPading (SNR = 20GB 22 aso ceceseesal and pelles haste naventeencaeeae ssn sutees agangass 54 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat 
Ricean Fading (SINK = Dvd B)ccsc2 thie cates gots eded ss vances cs etcecececeaiusl secdeash a euitiatee 54 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Selective Rayleigh Fading (SNR = 20 dB). .......eeceeeeeeeesseceesteeeeneeeeeeeeeeneeeenas aD 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Selective Rayleigh Fading (SNR = 5 dB). .........eeceeseeecceeeeeseetenceeneesoneesseteenee 55 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Selective Ricean Fading (SNR = 20 dB). oo... eee ceeeeecesececeseeeceeeeeeseeeeneerenes 56 
Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency- 
Selective Ricean Fading: (SNR = SUB). .4,.ccedasedoseeceshccgeeshceciseos iansepasnctavseneaeens 56 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency- 
Selective Rayleigh Fading (SNR = 20 CB)... eeeeeeeseeeseeeeeeeneecnaeeneeeeees a7 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency- 
Selective Rayleigh Fading (SNR = 5 dB)... cee eeeeeeeeeeeceseeeeeeeeseecnneenseeesees 57 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency- 
Selective Ricean Fading (SNR = 20 CB). ......ee ce eeeeceeceeesececeseeeceeeeeseeeeeneeenes 58 
Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency- 
Selective Ricean Fading (SNR = 5 AB). ......eeeeeecceeseceeceeececseeeecseceeeseeeeeneeeesaes 58 
Cumulants Expressed In Terms of Moments. .............:cescesseeeeseecsseceseeeeeeeeseenes 61 


XX 


EXECUTIVE SUMMARY 


Blind modulation classification remains a challenging problem despite the 
numerous studies that have investigated it. This thesis applies higher-order moments and 


cumulants as features to several modulation classification algorithms. 


An overview of digital communications is provided along with the statistical 
models for noisy and fading channels. The extraction of moments and cumulants is 
discussed at length, and the criteria for selecting them as features for discrimination are 
explained. Attempts are made to make the classifers robust to fading effects, first by 
investigating the effects of phase shifts on the statistics, and second by applying realistic 
power normalizations. Principal Component Analysis (PCA), Linear Discriminant 


Analysis (LDA) and their nonlinear, kernel-based counterparts are presented. 


Principal Component Analysis works slightly better than Linear Discriminant 
Analysis in most channel conditions and equals the performance of kernel-based PCA 
with less computational overhead. Results show that this classifier works very well with 
signals that are corrupted only by additive white Gaussian noise (AWGN). It achieves a 
classification rate of 95.1% at an SNR of 20 dB and 89.0% at 5 dB. In moderate fading 
conditions, such as might be experienced by a mobile device being carried by a person 
walking, it performs nearly as well, with classification rates in a Rayleigh channel of 
87.0% (20 dB) and 86.8% (5 dB), and 91.5% (20 dB) and 91.7% (5 dB) in a Ricean 
channel. Even PCA’s performance degrades sharply, however, in more severe fading 


conditions. 
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I. INTRODUCTION 


A. OBJECTIVES 


Blind modulation classification in which the modulation type of a waveform is 
determined without a prior! information has many potential applications. In the 
commercial sector, it would enable the interoperability of diverse communications 
devices. In the military domain, it would allow the demodulation of an intercepted signal 
to be performed in situ rather than after the signal has been sampled and forwarded on for 
processing. Software-defined radio would also benefit from this capability since the 
transmitter and receiver would not have to agree on a particular modulation scheme in 


advance. 


Because of the obvious benefits of blind modulation classification, much has been 
published on this topic in the literature in recent years, as will be seen below. There is no 
widespread agreement, however, on which modulation schemes to consider. Even within 
the main families of schemes, there are endless variants. Some papers consider phase- 
shift keying (PSK) versus frequency shift keying (FSK). Others introduce quadrature 
amplitude modulation (QAM) or variants of PSK such as minimum-shift keying. This 
thesis is no different in its arbitrariness. Nine different modulation schemes are 
considered: 2, 4, and 8-PSK, 2, 4, and 8-FSK, and 16, 64, and 256-QAM. How the 
features used for classification were chosen and the classifier algorithms are discussed. 
Simulation results for signals in additive white Gaussian noise (AWGN) and signals in 


fading channels are also discussed. 


B. BACKGROUND 


Various methods have been proposed to accomplish blind modulation 
classification. Marchand proposes the use of higher-ordered moments and cumulants in 
[1]. In one way or another, these features are used in [2, 3, 4]. Hatzichristos uses them as 


inputs to a neural network classifier [2], whereas Young uses a simpler approach 


involving threshold values [3]. Like, Chakravarthy, Ratazzi, and Wu propose a neural 
network-based classifier that uses cyclostationarity detection in its first tier and 4" and 


8"_ordered cumulants thereafter [4]. 


Other methods involve the use of wavelet transforms. Ho, Prokopiw, and Chan 
propose a scheme that uses very few symbols (on the order of 100) but do not consider 
QAM, in which both the amplitude and the phase are modulated [5, 6]. Prakasam and 
Madheswaran extracts several statistics from the wavelet transform and places them in an 


elaborate decision tree [7]. 


A summary of many of these methods can be found in [8], in which four broad 
categories for all the methods to date are defined: algorithms based on instantaneous 
amplitude, phase, and frequency; algorithms based on the wavelet transform; algorithms 


based on cumulants; and algorithms based on cyclostationarity properties. 


This thesis falls into the third category. Some of the same features — higher- 
ordered moments and cumulants - investigated in [2, 3] are examined, but a different set 
of classification algorithms is investigated. Specifically, the classification rates of 
Principal Component Analysis (PCA) and Linear Discriminant Analysis (LDA), as well 
as their nonlinear (kernel-based) counterparts, are determined and compared. These 
schemes are well-known and have been applied successfully to face-recognition 


applications in [9, 10]. 
C. ORGANIZATION 


The theory of digital communications systems, the nine modulation schemes 
considered in this work, and the mathematical models for introducing noise to wireless 
channels are reviewed in Chapter Hl. The features used for classification, namely 
moments and cumulants, and the reasons for choosing the particular moments and 
cumulants used in this work, are discussed in Chapter II]. The linear and nonlinear 
classification schemes that use the moments and cumulants as inputs are introduced in 
Chapter IV. The software implementation of the classifiers and the simulation results are 
discussed in Chapter V. Finally, conclusions and suggestions for further work are 


presented in Chapter VI. 


D. SIMULATION SOFTWARE 


All simulations were conducted in MATLAB version 7.6. Various functions from 
the Communications Toolbox were used. The code developed for this thesis is listed in 
the appendices. Note that this work implicitly assumes that some preprocessing has 
already occurred to down-convert the received signal as the functions in the 


Communications Toolbox used in this work simulate digital signals at baseband. 
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Yl. DIGITAL MODULATION AND WIRELESS CHANNELS 


The reader is introduced to the digital modulation schemes considered in this 
thesis in this chapter. Additionally, the mathematical models used to represent modulated 


signals corrupted during transmission are presented. 
A. THE NEED FOR MODULATION 


The purpose of any digital communication system is to transmit data from an 
information source to an information sink. At a minimum, a baseband waveform must be 
constructed from the symbols making up the information. This often amounts to 
assigning two different voltage levels to the binary symbols that represent the 


information. Construction of this baseband waveform is called pulse code modulation. 


Some media, such as coaxial cable, can propagate baseband waveforms with no 
further manipulation required. For wireless transmission, on the other hand, baseband 
waveforms are impractical for several reasons [11]. First, the antenna size necessary to 
transmit an electromagnetic wave is inversely proportional to the wave’s frequency. 
Therefore, shifting the baseband waveform to a higher frequency allows for a smaller 
antenna. For many mobile devices, this results in a carrier frequency on the order of 1 


GHz. 


Second, multiple transmitters operating at the same frequency generally interfere 
with each other. By assigning each transmitter a unique portion of the electromagnetic 
spectrum, multiple transmitters can operate without interference. This is called 


frequency-division multiple access. 


The translation of a baseband waveform to a higher frequency is accomplished by 
manipulating the features of a sinusoidal wave at that frequency. This process is called 


bandpass modulation. 


B. DIGITAL MODULATION SCHEMES 


1. Introduction 


There are three features of a sinusoidal wave that the designer can manipulate: 
amplitude, frequency, and phase. A given modulation scheme works on one or more of 


these. This work will consider three fundamental schemes: PSK, FSK and QAM. 
2. Phase Shift Keying 


As the name implies, PSK works by assigning a unique phase to each symbol 
present in the baseband waveform. Thus, binary phase-shift keying (BPSK) in which 
each symbol represents one bit uses two different phases, 4-PSK uses four different 


phases, and so on. The resulting waveforms can be represented mathematically by [12]: 





2E 20 
s(t)= S cos| 27 ft+—i| i=0,1,....M—-1, 2.1 
:(t) | 7 f: 7 (2.1) 


S 


where M is the number of bits per symbol, f. is the carrier frequency, E, is the energy 


per symbol, and 7) is the symbol interval. 


The case where M =4 illustrates a result that is widely exploited in 


communications systems: orthogonal signaling. The four phases are 0,%4,7, and 3%. 


Note, however, that 
a : 32 : 
cos [aes + =) =sin(27f,t) and cos [aes + =) =-sin(2af,t). 


This is important because the sine and cosine waves are orthogonal signals, 
meaning that they are uncorrelated in time over a symbol duration. Even when their 
frequencies are the same, they can be transmitted by the same antenna and recovered by 
the same receiver without mutual interference. When two signals are orthogonal, it is 
common to represent them as two perpendicular vectors that form a basis in R*. This 
leads to a conceptualization of the modulated symbols as points in a constellation, as 


shown in Figure 1. 





Quadrature Channel 





In-Phase Channel 














Figure 1. QPSK Constellation. 


Thus, any modulated waveform can be constructed as the sum of two sinusoids. 
The first, represented by the horizontal axis, is called the in-phase (I) channel. The 
second, represented by the vertical axis, is exactly out of phase with the first and is called 
the quadrature (Q) channel. For this reason, 4-PSK is usually referred to as quadrature 
phase-shift keying (QPSK). This result will be particularly important in the discussion of 
QAM. 


PSK schemes are popular for their symbol error performance and bandwidth 
efficiency. This comes at the price, however, of increased receiver complexity. In order 
to recover the baseband signal, the receiver must generate a sinusoid that is phase- 
synchronized with the received signal (except for one variant called differential phase- 
shift keying, which is not considered in this work). Modulation schemes that impose this 


requirement on the receiver are said to be coherently detected. 


PSK signals sometimes have an additional phase offset that rotates the 


constellation by a uniform amount, typically 74 for BPSK and 7% for QPSK. As will 


be seen in Chapter III, this can affect the features used for classification. 
3. Quadrature Amplitude Modulation 


Whereas PSK only uses phase differences to distinguish between symbols, QAM 
modulates both the phase and the amplitude of the carrier. Thus, it can be thought of as a 
generalization of PSK where all symbols need not have equal energy. Like PSK, the 
mathematical expression for an M-ary QAM signal can be decomposed into orthogonal 


sinusoids [12]: 


s,(t)= Pea, cos(2nf.t)- 2p sina ft) i=ObnnM=1 22) 


Figure 2 shows the symbol constellation obtained for a 16-QAM signal. 
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Figure 2. 16-QAM Constellation. 
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Although this work only considers QAM schemes that have square-shaped 
constellations (i.e., an equal number of symbols occupying each quadrant), this is not a 
general requirement for QAM. QAM signals tend to be bandwidth efficient and, like 
PSK, must be coherently detected. 


4. Frequency-Shift Keying 


FSK signals are unlike PSK and QAM in that each of the M symbols in an M-ary 


scheme has its own signaling frequency [12]: 





s,(t) = - cos(27 ft) i=1,2,...,M. (2.3) 


FSK uses bandwidth inefficiently but can be detected noncoherently, which 
reduces receiver complexity. It cannot be decomposed and visualized as I- and Q- 


channels, since it uses multiple signaling frequencies. 


C. SOURCES OF SIGNAL DEGRADATION 


1. Introduction 


In the process of transmission and reception, there are several ways in which a 
signal can undergo degradation, making it harder for the receiver to demodulate it 
correctly—or making it harder to determine its modulation type. A radio frequency (RF) 
signal always experiences a decrease in power due to the distance from the transmitter, 
called path loss. In addition, this work considers two types of degradation: AWGN and 
multipath fading. 


2 Additive White Gaussian Noise 


Any receiver will have a certain amount of thermal energy associated with its 
antenna and other circuitry. Additionally, interfering signals will be received by the 
antenna [11]. This interference can be modeled as a Gaussian random process that adds 


to the modulated signal. For a signal propagating in free space, this model is sufficient to 


describe the degradation of the signal. It is common to express this process in terms of its 


. N os F 
power spectral density, 9 4 , which is a constant as a function of frequency. 





N/2 
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Figure 3. Power Spectral Density of AWGN. 


The effect of AWGN on a signal is to introduce uncertainty as to which symbols 
are being transmitted. The variance in the thermal noise will shift each symbol around its 
intended value. For example, the constellation of a 16-QAM signal with 10,000 symbols 
in the presence of AWGN is shown in Figure 4. 
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Figure 4. Effect of AWGN on 16-QAM (SNR = 20dB). 


This degradation can be quantified by defining the signal-to-noise ratio (SNR), 
which is simply the ratio of the signal power to the noise power. For the purpose of error 


analysis, it is necessary to modify this quantity to account for symbol time and receiver 


bandwidth. This results in the quantity “yy , which denotes the ratio of symbol energy- 
0 
to-noise energy. Note that much of the literature in communications theory expresses 
SNR as % , the ratio of bit energy to noise energy. This is because the usual measure 
0 
of error performance is bit error ratio, and expressing it in terms of Be allows for an 


apples-to-apples comparison between M-ary systems. In order to keep the bit energy 
consistent, however, an M-ary system must have M times as much symbol energy as a 


binary system (E,=ME,), and, therefore, M times as much power, assuming equal 
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symbol rates. Because the goal of this work is to distinguish between modulation types, 


the symbol energy is kept consistent and oA is used as the measure of SNR. 
0 


3. Fading Channels 


While AWGN alone adequately models the degradation of a signal in free space, 
many signals undergo more complicated propagation on top of AWGN. The presence of 
physical obstructions such as buildings causes multipath propagation losses. As the name 
implies, this means that the signal arrives at the receiver via more than one propagation 
path. Depending on the relative lengths of the paths, this can result in either constructive 
or destructive interference. Furthermore, relative motion between the transmitter and 
receiver causes Doppler shift. 

Multipath propagation results in both large scale and small scale fading. Large 
scale fading is similar to the path loss that a signal experiences in free space in that it is 
proportional to Vn , where r is the distance between the transmitter and receiver. In free 
space n=2, whereas in a fading channel, n depends on the particular objects reflecting 
the signal. Also, this relationship only determines the mean path loss due to large scale 
fading. The actual path loss can be modeled as a random variable with a log-normal 
variation around this mean value [11]. 

Small scale fading refers to finer variations in the relative position of the 
transmitter and receiver that can cause multipath copies of the signal to interfere with 
each other. When a line-of-sight component is present among the reflected paths, the 
received amplitude can be modeled as a random variable with a Ricean probability 


density function (PDF) [11], 
I, r+A’ rA 
P(n)=— e) 552 Is ™ } (2.4) 


where 7, is the distance between transmitter and receiver, A is the magnitude of the line- 








of-sight component, o is the mean amplitude of the reflected components, and /, is a 


modified Bessel function of the first kind. When no line-of-sight component is available, 


this PDF simplifies to the Rayleigh PDF: 
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I, r 
o(a)=ieexe|~ >} (2.5) 


Finally, in mobile applications the Doppler effect causes the multiple received 
copies of the signal to alternate between constructive and destructive interference. When 
this effect is slow compared to a symbol interval, it is called slow fading. Similarly, 


when the effect occurs within a symbol interval, it is called fast fading [13]. 


Fortunately, the MATLAB Communications Toolbox includes functions that 
account for all these effects. A 16-QAM signal in slow Rayleigh fading is shown in 


Figure 5. The overall effect is to rotate the constellation. 
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Figure 5. Effect of Rayleigh Fading on 16-QAM (SNR = 20dB). 


D. CONCLUSION 


In this chapter, the reader was introduced to some of the digital signaling 
techniques commonly in use and the statistical methods that model their corruption by 
white noise and multipath fading, which both make the task of classifying the signals 


harder. The features used as a basis for discrimination will be discussed in Chapter III. 
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Hl. SIGNAL FEATURES 


In order to implement the classifier, features must be extracted from the signals to 
provide a basis for distinguishing between the modulation types. In this chapter, the 
features that will be used in classification are presented and their behavior under AWGN 


and multipath fading examined. 
A. INTRODUCTION 


For FSK signals, the most effective approach is to exploit their unique spectral 
characteristics. For PSK and QAM, a number of different discriminating features have 
been proposed in the literature. The use of higher ordered moments and cumulants was 
proposed in [14] and explored in [2, 3]. Wavelet transforms were examined in [5, 6, 7]. 


More recently, [4] presented an approach based on cyclostationarity concepts. 
B. SPECTRUM BASED CLASSIFICATION OF FSK SCHEMES 


The spectrum of an FSK signal is sufficient to distinguish it from PSK and QAM 
signals and to separate different M-ary FSK schemes. A simple way to automate this 
examination of the spectrum is presented in [15]. First, the fast Fourier transform (FFT) 
of the signal is computed. Next, transform peak values are compared to each other in 


order to determine whether the signal is FSK and if so, what type. 


Consider the single-sided spectrum of the BFSK signal shown in Figure 6. There 
are clearly two peaks in this spectrum. If the ratio of the second highest value to the third 
highest value in the spectrum is computed, the result is a number larger than one. 
Likewise, a 4-FSK signal will have a large ratio of the fourth to the fifth highest value 
and an 8-FSK signal of the eighth to the ninth highest value. For non-FSK signals, these 


ratios will all be approximately unity. 
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Figure 6. Single-Sided Spectrum of BFSK. 


C. MOMENTS 


Statistical moments are nothing more than the expected value of a random 
variable raised to the power indicated by the order of the moment. A first order moment 
is the mean of the variable. In many applications, the second order moment is a 
measurement of power. For complex random variables, a conjugate term is often 
included (note that the signals generated in this work are simulated by their complex 


envelope). In general, the moment of a random variable is defined as 
E,,5=E|x*(z)’|, a,beZ, 3.1) 


where x is the variable and x is its complex conjugate. The order of this moment is 
a+b. For our purposes, x is a vector of samples of a digitally modulated signal. Rather 
than explicitly calculating statistical moments, they are estimated by raising each element 


in the vector to the indicated power and computing the numerical mean. In this work the 
16 


moments are normalized by the measure of signal power appropriate for the order of the 
moment. For example, a second order moment is normalized by the signal power itself, a 
fourth order moment is normalized by the power squared, and so on. Additionally, all 


moments are central moments as the signal mean is zero. 
D. CUMULANTS 


In order to define the cumulants of a random variable x, we first define its 


characteristic function and second characteristic function, respectively, as [16]: 


—00 


O(s)= | f(x)e"dx (3.2) 


and 
Y(s)=n@®(s). (3.3) 
The n'-order cumulant of x is defined as the n™ derivative of the second 


characteristic function evaluated at s = 0: 


422 (3.4) 
ds" 


Note, that x may be complex-valued. In order to be consistent with the notation 





above for moments, cumulants will be represented by C 


x,a,b°* 


Fortunately, cumulants can also be expressed as functions of equal and lower 
ordered moments, making it simple to compute them numerically. These expressions 


were computed up to the eighth order in [1] and are presented in Appendix A. 
E. FEATURE SELECTION 


Some points should be made about the application of moments and cumulants to 
the problem at hand. First, they will have imaginary components. It was asserted in [1] 
and [2] that the imaginary components would be zero for the signals considered here. 
This may be the case if there is no rotation of the signal constellation. As mentioned in 
Chapter II, however, fading channels have the effect of rotating the entire constellation; 


also, some implementations of BPSK and QPSK may include an arbitrary phase offset. 
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To illustrate the effect of a phase rotation, E,,, is calculated below for a QPSK 


signal with an arbitrary phase shift of ¢. The expression used is the baseband (or 


complex-envelope) expression for QPSK. Since QPSK signals have four possible values, 


d(k) is a uniformly distributed random integer between 0 and 3. 


e [oo{ aeres))| |-e[oo[ (Sac) 


eal 
= exp(j4g) > 7 exp j2ak) (3.5) 


k=0 
= exp( jag) [1414141] 


= exp(j4¢). 

Clearly, the phase rotation may introduce an imaginary component into the 
Statistic that varies with the phase shift ¢. Section A of Appendix C contains a general 
proof that, for PSK and QAM signals, an arbitrary phase shift of the constellation affects 
the phase of the moments but not their magnitudes. In section B, this result is shown to 


hold for the cumulants C,,9,C.11> Cxao» Ce3i> C602 Crs: and C,g . Therefore, the 


magnitude of each of these statistics are used rather than the real component alone. In 
Tables 1 and 2 the magnitudes of the moments and cumulants computed for all nine 
modulation types using uncorrupted signals of 20,000 symbols each are presented. These 
were estimated by simulating the signals in MATLAB. Note that they have been 


normalized by the estimated received signal power. 
The main purpose of presenting the moments and cumulants here is to verify that 


they are being calculated correctly, but it should already be clear that some of them will 


be more useful than others. For example, C,,, is useful for identifying BPSK, but 


confuses all the other modulation types. In contrast, C,,, and C,;, have different values 


for all PSK and QAM schemes, although the 64- and 256-QAM values are very close 
together. As will be seen, adjustments are necessary in order to make them more robust 
to noisy conditions. Accordingly, their behavior in noisy conditions is investigated 


before deciding which ones to use as classification features. 
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Table 1. 


Moments of Digitally Modulated Signals (SNR=~ ). 
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BPSK | QPSK | 8PSK | BFSK | 4FSK | 8FSK | 16QAM | 64QAM | 256QAM 

Foo 
1.000 0.006 | 0.006 | 0.000 | 0.000 | 0.000 0.008 0.007 0.008 

Fai 
1.000 1.000 | 1.000; 1.000; 1.000 | 1.000 1.000 1.000 1.000 

Fao 
1.000 1.000 | 0.006 | 0.000; 0.000 | 0.000 0.681 0.618 0.606 

| ee 
1.000 0.006 | 0.006 | 0.000 | 0.000 | 0.000 0.011 0.012 0.013 

Fos 
1.000 1.000 | 1.000; 1.000; 1.000 | 1.000 1.320 1.381 1.396 

Eg 
1.000 0.006 | 0.006 | 0.000 | 0.000 | 0.000 0.018 0.025 0.025 

E54 
1.000 1.000 | 0.006 | 0.000 | 0.000 | 0.000 1.321 1.296 1.291 

Fao 
1.000 0.006 | 0.006 | 0.000 | 0.000 | 0.000 0.019 0.022 0.025 

E33 
1.000 1.000 | 1.000 | 1.000 | 1.000 | 1.000 1.960 2.225 2.293 

E80 
1.001 1.000 | 1.000 | 0.000; 0.000 | 0.000 2.204 1.907 1.828 

Ea 
1.001 0.006 | 0.006 | 0.000 | 0.000 | 0.000 0.031 0.050 0.054 

E62 
1.001 1.000 | 0.006 | 0.000 | 0.000 | 0.000 2.485 2.755 2.815 

F543 
1.001 0.006 | 0.006 | 0.000 | 0.000 | 0.000 0.032 0.047 0.055 

Fas 
1.001 1.001 | 1.000 1.000 | 1.000 | 1.000 3.124 3.961 4.194 





Table 2. | Cumulants of Digitally Modulated Signals (SNR=°). 



























































BPSK | QPSK | 8PSK | BFSK | 4FSK | 8FSK | 16QAM | 64QAM | 256QAM 

CL 20 
1.000 0.006 | 0.006 |} 0.000 | 0.000 | 0.000 0.008 0.007 0.008 

Co44 
1.000 1.000 | 1.000 | 1.000} 1.000 | 1.000 1.000 1.000 1.000 

Co40 
2.000 1.000 | 0.006 | 0.000} 0.000 | 0.000 0.681 0.618 0.606 

Co31 
2.000 0.012 | 0.012 | 0.000 | 0.000 | 0.000 0.012 0.011 0.011 

CL. 2. 
2.000 1.000 | 1.000 | 1.000} 1.000 | 1.000 0.680 0.619 0.604 

C.60 
15.999 0.087 | 0.006 | 0.000 | 0.000 | 0.000 0.068 0.054 0.052 

Co 51 
15.999 3.998 | 0.024 | 0.000} 0.000 | 0.000 2.084 1.794 1.738 

CLs 2 
15.999 0.062 | 0.065 | 0.000 | 0.000 | 0.000 0.060 0.050 0.048 

C33 
15.999 4.000 | 4.000 | 4.000} 4.000 | 4.000 2.080 1.798 1.733 

C,.80 
243.962 | 33.971 | 0.998 | 0.000} 0.000 | 0.000 14.023 11.465 11.009 

Coa 
243.957 1.087 | 0.012 | 0.000 | 0.000 | 0.000 0.870 0.715 0.737 

Cy 2 
243.962 | 45.984 | 0.280 | 0.000} 0.000 | 0.000 29.864 27.043 26.478 

C53 
243.962 0.396 | 0.528 | 0.000 | 0.000 | 0.000 0.273 0.328 0.345 

CoA 
243.962 | 17.999 | 16.992 | 17.000 | 17.000 | 17.000 17.371 24.098 25.728 


One potential problem in using these statistics as they are is that the magnitude of 
the cumulants increases with their order. This characteristic could have the unintended 


consequence of weighting these larger statistics more heavily in the classification 
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scheme. 


In order to mitigate this effect, [17] proposed raising each cumulant to the 


power 2%, where n is the cumulant’s order. 


cumulants is presented in Table 3. 


A revised table of these normalized 









































Table 3. Normalized Cumulants of Digitally Modulated Signals (SNR=00). 

BPSK | QPSK | 8PSK | BFSK | 4FSK | 8FSK | 16QAM | 64QAM | 256QAM 
CL20 

1.000 | 0.007 | 0.006 | 0.000 | 0.000 | 0.000 0.007 0.007 0.008 
Co44 

1.000 1.000 | 1.000 | 1.000; 1.000 | 1.000 1.000 1.000 1.000 
Co 40 

1.414 | 1.000 | 0.070 | 0.000 | 0.000 | 0.000 0.824 0.787 0.778 
Co31 

1.414 | 0.107 | 0.103 | 0.000 | 0.000 | 0.000 0.098 0.100 0.099 
Co 20 

1.414 1.000 | 1.000 | 1.000; 1.000 | 1.000 0.824 0.786 0.777 
C..6.0 

2.520 | 0.424 | 0.173 | 0.004 | 0.004 | 0.004 0.368 0.369 0.372 
C51 

2.520 1.587 | 0.265 | 0.006 | 0.006 | 0.006 1.276 1.216 1.201 
Co 42 

2.520 | 0.379 | 0.384 | 0.006 | 0.006 | 0.007 0.351 0.356 0.350 
C33 

2.520 1.587 | 1.587 | 1.587 | 1.587 | 1.587 1.276 1.215 1.200 
C,80 

3.952 2.414 | 0.999 | 0.018 | 0.018 | 0.018 1.933 1.842 1.820 
Co 

3.952 | 0.982 | 0.313 | 0.018 | 0.018 | 0.018 0.895 0.897 0.905 
Co 62 

3.952 | 2.604 | 0.676} 0.040 | 0.039} 0.040 2.337 2.282 2.268 
C53 

3.952 | 0.763 | 0.818 | 0.036 | 0.036 | 0.037 0.693 0.719 0.768 
CoA 

3.952 2.060 | 2.030 | 2.031 | 2.031 | 2.031 2.042 2.218 2.254 
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Simulations showed that taking the magnitude of the statistics and normalizing 
the cumulants according to their order greatly improved discrimination power of the 
features considered when dealing with AWGN. However, simulations also showed that 
some of the statistics were very sensitive to the received signal power. In addition, the 
fading process affects the power of signals as well as their phase. A simple assumption 
will improve the classifier’s performance in AWGN and allow it to deal better with faded 
signals. Suppose that a measurement of received power is available in the absence of any 
signals. This measurement will equal the power added to the signals by AWGN. 
Subtracting this from the received signal power will allow us to normalize the statistics 


by the “noise-free” power. 


As stated previously, all FSK signals have nearly the same moments and 
cumulants and are handled by measuring their spectral characteristics. For the remaining 
six modulation types, the number of features clearly must be reduced in order to simplify 
the problem for the classifier. Thus, the challenge is to identify a smaller set of moments 
and/or cumulants that can distinguish between the modulation schemes and that are 


relatively invariant over a range of SNRs. 


In order to accomplish this task, MATLAB was used to generate and corrupt the 
modulated signals with AWGN and then to extract their moments and cumulants for even 
orders from two to eight. These results are plotted for each statistic in Section A of 
Appendix B. Sections B through I show the same statistics under Rayleigh and Ricean 
fading conditions for all combinations of slow or fast and frequency-flat or frequency- 


selective channels. The exact parameters used in MATLAB are also presented. 


Sections B and C show that the moments and cumulants do not change much 
under slow, frequency-flat fading conditions for either Rayleigh or Ricean channels, 
indicating that the classifier should work well in these regimes. In sections D and E, the 
frequency-selective nature of the channels begins to shift the values, although some of 


them are still stable with increasing AWGN. 


In sections F and G, the values have shifted significantly due to fast fading. 


Furthermore, they are less stable with increasing AWGN. Since these results are for 


Ze, 


frequency-flat channels, we can conclude that Doppler shift has a more significant effect 
on the statistics than path delay. In sections H and I, which show the effects of fast, 
frequency-selective fading, the statistics are somewhat more stable but still well off their 


original values. 


Three statistics were selected as classification features: C,,,,C,;,,and C,,). All 


x,5,1? 
of them are invariant to phase shifts and provide good separation between the PSK 
schemes and some separation between the QAM schemes. Their plots are shown in 


Figures 7 to 9. While some of the moments (such as E,,, and E,,,) appeared to meet 


the same criteria, including them tended to make the classifier perform worse with faded 
signals. The values for 64- and 256-QAM do not differ until the second or third decimal 


place, indicating that they will present a challenge for the classifier. 
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Figure 7. C..4.9 (AWGN only). 
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C,;, (AWGN only). 
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Figure 9. C..¢9 (AWGN only). 


Figures 10 through 12 show the same cumulants for a slow, flat Rayleigh channel 


while Figures 13 through 15 show them for a fast, frequency selective channel. The 


values do not change much under moderate fading conditions, but under more severe 


fading they change significantly. Since any classifier is only as good as the features it 


works with, this indicates that severe fading conditions will be problematic. 
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Figure 10. 


C.49 (slow, flat Rayleigh fading). 
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Figure 11. 


C,;, (slow, flat Rayleigh fading). 
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Figure 12. 


C29 (slow, flat Rayleigh fading). 
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Figure 13. 


C..49 (fast, frequency-selective Rayleigh fading). 
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Figure 14. 


C,,;, (fast, frequency-selective Rayleigh fading). 
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Figure 15. 


C..¢9 (fast, frequency-selective Rayleigh fading). 





F. CONCLUSION 


So far, we have seen the features that will form the basis for classification. These 
features will be grouped into 3-dimensional vectors for exploitation by different 


classifiers, which are the subject of Chapter IV. 
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IV. LINEAR AND NONLINEAR CLASSIFIERS 


Now that discriminating features have been selected, a classification scheme must 


be identified to exploit them. This chapter discusses several classification methods. 
A. INTRODUCTION 


The use of higher-ordered moments and cumulants to classify modulation types 
was explored in [2] and [3] using a neural network and a classification tree, respectively. 
The purpose of this thesis is to investigate a different set of classification schemes, all of 
which rely on eigenvector decomposition in some way. This chapter will present two 


linear methods, PCA and LDA, and their nonlinear counterparts. 
B. LINEAR CLASSIFICATION SCHEMES 


Linear classification schemes work by projecting data into a feature space using a 
linear mapping and then comparing the result to a centroid for each class. If the data is 
linearly separable, these schemes work well. These algorithms attempt to find the best 


projection matrix. 
1, Principal Component Analysis 


PCA seeks the best representation of the data in a least-squares sense [18]. It does 
this by decomposing the data covariance matrix into its eigenvectors and choosing the 
most significant of them to form a projection matrix. The training data in this case is a 
matrix X of column vectors, denoted x,, containing the chosen signal statistics. Assume 
that there are a total of P of these statistical profiles. Let N be the number of classes (i.e., 
modulation types) represented in the training set. Also let n be the number of training 


profiles in each class so that P=nN. 


Define the mean profile as 


m=+> x, (4.1) 


The first step is to subtract the mean profile from each of the P profiles since it is 


common to all of them and thus contains no information useful for classification: 


x, =X,—m. (4.2) 


L 


These new centered training profiles form the new centered training matrix x. 
Define the covariance matrix as 
C=xXXx' (4.3) 
and decompose it into its eigenvectors w,. The eigenvectors corresponding to the top k 
eigenvalues form the projection matrix W. In this application, all of the possible 
eigenvectors were significant and were included in the projection matrix. 
The columns of W define the feature space into which the statistical profiles will 
be projected. In order to determine the centroids for each class, first project the centered 


training matrix into this space: 
T =W'X =|W' x IW? x21... WT xe |. (4.4) 


The centroid for class / will be 
1< = 
2, = a (4.5) 
i=] 


In order to use the classifier, compute the projection of a test profile ¢ into the 
feature space, 
t'=W' (t—m), (4.6) 
and then choose the centroid having the smallest Euclidean distance to the point r'. 


2. Linear Discriminant Analysis 


In one sense, PCA is just a compression scheme in that it constructs the feature 
space by arbitrarily selecting the most significant eigenvectors as basis vectors. LDA, on 
the other hand, seeks a different projection matrix W that will maximize the separation 


between classes [18]. 


To explain this further, it is necessary to define the within-class scatter matrix 


S,, and the between-class scatter matrix S,, defined as 
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5, =>'S, (4.7) 


i=l 


where S, is the covariance matrix of class i, (defined in the same way as Equation 3.13 


above, except that it only contains profiles from class 7) and 
N 
Sy = din, —m)(m, —m)', (4.8) 
i=l 


where m, is the mean profile of class i and m is defined as in Equation 3.11. If these 
same matrices were calculated in the feature space formed by projecting the profiles 
using the projection matrix W, the results would be W’S,,W and W’S,W, respectively. 
The problem for LDA becomes finding the matrix W that maximizes the ratio of the 
determinants of these two matrices, 


T 
ew (4.9) 
From [18], finding the matrix W in the expression above is equivalent to solving 
the Eigen problem 
S,W, =A,SyW,;, (4.10) 
where w, is the ith column of W for i=1 to the number of features. As with PCA, the 
columns are arranged in order of decreasing eigenvalues and an arbitrary number are 
kept. 
Note, in some applications it is necessary to apply PCA before LDA; otherwise, 
S, may be singular. This requirement is typical in applications where the dimension of 
the data is larger than the number of training observations. However, the dimension of 
the feature vector is three in our application, as only three representative features were 


selected, and LDA can be applied directly. 


The class centroids and the projection of a given signal profile into the feature 


space are calculated in the same manner as in PCA. 
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C. NONLINEAR CLASSIFICATION SCHEMES 


PCA and LDA work well with classes that are linearly separable. To visualize 
linear separability, consider a projected feature space which has three dimensions. If the 
classes are linearly separable, simple planes can be drawn as the decision boundaries 
between the classes. For data that is not linearly separable, it would be preferable to draw 
curved planes separating the classes. Conceptually, this can be achieved by a nonlinear 
mapping of the data into a higher dimensional space in which the classes are linearly 
separable. In the original space, the decision boundaries between them would appear 


curved. 
1. The Kernel “Trick’”’ 


Carrying out an arbitrary nonlinear mapping to the higher-dimensional space may 
be computationally expensive. However, if the mapping function is chosen well, it may 
be encapsulated in a mechanism commonly known as the kernel “trick.” Rather than 
operating on individual data points, kernel functions compute the dot product between 
two data points after projection into the higher dimensional space. The projection of an 
individual data point into the space is formed by computing its dot product with every 


other member of the space and arranging the dot products in a vector. 


As an example, consider the following mapping function from two to three 


dimensions presented in [19]: 
(%,.)) > (27,53, V2x J. (4.11) 
The goal is to find a function that will express the dot product between two points 


x and y in the higher dimensional space in terms of a function in the lower dimensional 


space. It can be shown that the function 


k(x,y)=(x-y) (4.12) 


accomplishes this goal. Common kernel functions include the following: 
polynomial: 
d 
k(x, y)=(a(x-y)+b) (4.13) 
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Gaussian: 


-_ — 2 
k(x, y) = exp [=e=p) (4.14) 
20 
sigmoid: 
k(x, y) = tanh(x«(x- y)+ 8). (4.15) 
2. Nonlinear Component Analysis 


Nonlinear Component Analysis is the extension of the PCA method to higher 
dimensional spaces via the kernel trick; hence, it is sometimes known as kernel PCA 
(KPCA). The first step is to construct the Px P kernel matrix K by applying the kernel 
trick to each pair of points in 64 (this can usually be done with a matrix equation). Thus, 
K is the matrix of training data in the higher dimensional space. Once again, the 
projection matrix W is formed by using the eigen-decomposition of the matrix K. If 
dimension reduction is desired, only the most significant eigenvectors need to be kept. 
However, the matrix K must be reformed in order for its dimensions to match other 
matrices by applying the following steps: 

K=WIw', (4.16) 


where I is the diagonal matrix of eigenvalues. At that point, the columns of W should 


1 
be normalized by the magnitude of each column in the matrix (w’ kw)” «oh: 


projected training data is calculated by 
T=W'K (4.17) 
and the class centroids are calculated in the same manner as with the linear PCA scheme. 


To project a test profile ¢, first translate it to the higher dimensional space by 


jaa x (4.18) 
and then project it for classification, leading to 
t"=W't'. (4.19) 
as General Discriminant Analysis 


As was the case for their linear counterparts, GDA differs from KPCA by taking 


into account the scatter within classes. The kernel matrix K is formed as in the previous 
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section and decomposed into eigenvectors U and eigenvalues [’. Next, the block 


diagonal matrix B is formed from N square matrices. The dimension of these matrices is 


nxn and each coefficient in the matrix is equal to |y,. 


These results are used to find the eigenvectors (@ of the matrix U’BU. The 


projection matrix is found by computing 


W=Ur'B (4.20) 
1 
and normalizing by the matrix (ww)? [20]. The class centroids and projected 


testing data are calculated in the same manner as the previous section. 


D. CONCLUSION 


Several well-known classification schemes were presented in this chapter. In the 


next chapter, their application to the features selected in Chapter III are presented. 
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V. IMPLEMENTATION AND RESULTS 


The previous four chapters have laid the groundwork for classifying digital 
signals. In this chapter, the generation of the signals is explained, and the 
implementation of the classifier is discussed. Finally, we present simulation results. All 


MATLAB code developed for this thesis can be found in Appendix D. 


A. SIGNAL GENERATION AND CORRUPTION 


The MATLAB Communications Toolbox provides standard objects for 
generating digital signals. As noted in Chapter I, these objects simulate the signals at 
baseband. For PSK and QAM, this means that the output is the complex envelope of the 
signal with one sample per symbol. By definition, FSK signals cannot be represented as 
a single frequency, so they are generated with several samples per symbol, but the 


resulting spectrum is still centered at 0 Hz. 


In terms of a physical system, this means that the digital signal has already been 
downconverted from its carrier frequency to baseband. In the case of PSK and QAM, it 
also means that the symbol rate has been obtained. Both of these assumptions are 


reasonable given an accurate spectrum measurement. 


For fading channels, the Communications Toolbox also provides functions to 
generate filter objects for Rayleigh and Ricean channels. For either type of fading, the 
number of paths must be specified along with the gain and time delay for each one. 
Additionally, the Doppler shift must be specified. These parameters were varied to 
simulate slow and fast fading as well as frequency-flat and frequency-selective channels. 
Whether or not the channel involves fading, AWGN can be added to the signals using a 
standard function in MATLAB. 


B. CLASSIFIER IMPLEMENTATION 


Separate functions were developed to implement the four classifiers discussed in 


Chapter IV, but they all peform similar tasks. First, fifty training signals of 20,000 


35 


symbols are generated for each modulation type and estimates for the three selected 
cumulants computed for each one. Next, these parameters are passed to a training 


function that determines a projection matrix and class centroids. 


In each trial of the testing phase, one thousand signals (also of 20,000 symbols) 
are generated and a modulation scheme is chosen at random for each one. One hundred 
trials are run for each combination of testing SNR and fading parameters in order to 
ensure that the results are statistically significant. The spectrum of each signal is 
analyzed to see if it matches one of the FSK types. To do this, an FFT is calculated on 
the first 4096 samples of the signal. The ratios of the second to the third peaks, the fourth 
to the fifth peaks, and the eighth to the ninth peaks are calculated and compared to 
threshold values estimated in a separate simulation. If the ratio of the second to the third 
peak is above the respective threshold, the signal is classified as BFSK. If not, similar 


comparisons are made for the other ratios to determine if the signal is 4-FSK or 8-FSK. 


If none of the thresholds is crossed, the signal is not FSK. Its cumulants are 
calculated and projected into the new feature space using the projection matrix calculated 
by the training algorithm. The projected features are compared to the class centroids and 
the closest one is chosen as the modulation type. A confusion matrix is saved so that the 
classifier’s performance can be analyzed for each modulation type. The procedure for 


classifying a test signal is shown in Figure 16. 
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Figure 16. Classification Flow Chart for Test Signals 


One question that must be addressed is the SNR at which the classifier is trained. 
A classifier trained at 20 dB will have better accuracy handling test signals at 20 dB than 
5 dB and vice versa. However, after applying the noise-free power normalization 
discussed in Chapter III, these differences were so slight that it was found to be simpler 


to train the classifiers with uncorrupted data. 
Ce LINEAR VERSUS NONLINEAR CLASSIFTERS 


The performance of the nonlinear classifiers was generally disappointing. As 
shown in Tables 4 through 11, for some kernels KPCA and GDA equaled the 
performance PCA and LDA, respectively. Note that in the case of KPCA, the best- 


performing scheme used a polynomial kernel of order one, which is just a more 
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Actual Modulation Type 


complicated implementation of PCA. Given that nonlinear schemes involve significantly 


more computations, the linear schemes should be used. 


The results that compare all four schemes are presented below as confusion 
matrices for signals in AWGN (more extensive results for the performance of PCA and 
LDA in fading channels is presented in section D). The overall performance and 
standard deviation are presented at the top of each table. Several different kernels were 
tried as well as different kernel parameters for each one. Only the kernel type leading to 
the best results is presented for KPCA and GDA schemes. PCA achieved the best overall 
classification rate, although section D shows that LDA performs better in fast Ricean 


channels. 
1. PCA 


PCA achieves the best overall performance and has the added benefit of being the 
simplest classifier to implement. At an SNR of 20 dB (shown in Table 4), the errors are 
primarily due to 64-QAM and 256-QAM being confused for each other. At 5 dB (shown 
in Table 5) there is more confusion between these two. Additionally, 16-QAM is 
sometimes mistaken for the other QAM schemes, and a significant number of 8-FSK 


signals are classified as 8-PSK 


Table 4. Confusion Matrix for PCA classifier in AWGN, SNR = 20 dB. 


Classifier Output (Average Performance = 95.1%, Standard Deviation = 0.6% 



































BFSK| QFSK| 8FSK| BPSK) QPSK} 8PSK) 16QAM| 64QAM/256QAM| Total 
BFSK 11197 0 0 0 0 0 0 0 0|11197 
QFSK 0} 11272 0 0 0 0 0 0 011272 
8FSK 0 0} 11016 0 0 106 0 0 0/11122 
BPSK 0 0 0} 11186 0 0 0 0 0|11186 
QPSK 0 0 0 0} 10903 0 0 0 0|10903 
8PSK 0 0 0 0 0} 10995 0 0 0|10995 
16QAM 0 0 0 0 0 0 11220 0 0|11220 
64QAM 0 0 0 0 0 0 5 9342 164010987 
256QAM 0 0 0 0 0 0 0 3166} —7952/11118 
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Actual Modulation Type 


Table 5. 


Confusion Matrix for PCA classifier in AWGN, SNR = 5 dB. 


Classifier Output (Average Performance = 89.0%, Standard Deviation = 1.1% 

































































QAM somewhat more and 256-QAM somewhat less. 


BFSK| QFSK| 8FSK| BPSK} QPSK) 8PSK/ 16QAM) 64QAMI256QAM| Total 

&|BFSK 11053 0 0 0 0 0 0 0 0|11053 
F |QFSK o| 10971 0 0 0 0 0 0 0|10971 
S|8FSK 0 1| 10865 0 0 452 0 0 0/11318 
= \BPSK 0 0 O| 11243 0 0 0 0 0/11243 
8 IQPSK 0 0 0 0} 11137 0 0 0 0/11137 
= 8PSK 0 0 0 0 0| 10967 0 0 0|10967 
2/16QAM 0 0 0 0 0 0| 10394 792 14]11200 
<64QAM 0 0 0 0 0 0 1371 5826] 3774|10971 
256QAM 0 0 0 0 0 0 489 4154) 6497111140 

2. LDA 


LDA performs as well as PCA at 20 dB (Table 6), although it misclassifies 64- 


perform as well as PCA due to misclassifying the QAM schemes more often. 


Table 6. 


Confusion Matrix for LDA classifier in AWGN, SNR = 20 dB. 


Classifier Output (Average Performance = 94.9%, Standard Deviation = 0.66%) 


At 5 dB (Table 7) it does not 















































BFSK| QFSK| 8FSK| BPSK) QPSK) 8PSK)  16QAM)  64QAM/256QAM| Total 
BFSK 11073 0 0 0 0 0 0 0 0|11073 
QFSK 0} 11102 0 0 0 0 0 0 0/11102 
8FSK 0 2| 11006 0 0 109 0 0 1611133 
BPSK 0 0 1} 11117 0 0 0 0 0/11118 
QPSK 0 0 0 0} 11309 0 0 0 0|11309 
8PSK 0 0 0 0 OQ} 11111 0 0 O11111 
16QAM 0 0 0 0 0 0 10946 0 0|10946 
64QAM 0 0 0 0 0 0 1 8220} 2892111113 
256QAM 0 0 0 0 0 0 0 2060) —_9035}11095 
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Table 7. 


Confusion Matrix for LDA classifier in AWGN, SNR = 5 dB. 


Classifier Output (Average Performance = 87.2%, Standard Deviation = 1.1%) 






















































































& BFSK| QFSK|  8FSK| BPSK) QPSK} 8PSK| 16QAM/| 64QAM/256QAM| Total 
Zz BFSK 11042 0 0 0 0 0 0 0 0111042 
L|QFSK 0} 10999 0 0 0 0 0 0 010999 
= 8FSK 0 1; 10774 0 0 432 0 0 2/11209 
6 (BPSK 0 0 0} 10984 0 0 0 0 0/10984 
= QPSK 0 0 0 0} 11145 0 0 0 0111145 
2 8PSK 0 0 0 0 0} 10955 0 0 5/10960 
<= /16QAM 0 0 0 0 0 0 9938 1226 9111255 
64QAM 0 0 0 0 0 0 1771 4847| 4538|11156 
256QAM 0 0 0 0 0 0 933 3829|  6488/11250 
3. KPCA 
For KPCA, the best kernel turns out to be k (x, y) = (x: y +1)! , which is actually 
a linear scheme. As one would expect, its performance is the same as that of PCA at both 
20 dB and 5 dB (Tables 8 and 9, respectively). Because of the extra computational 
overhead required by KPCA, PCA is preferred. 
Table 8. Confusion matrix for KPCA classifier in AWGN, SNR = 20 dB. 
Classifier Output (Average Performance = 95.1%, Standard Deviation = 0.7% 
BFSK| QFSK}  8FSK| BPSK} QPSK} 8PSK| 16QAM/| 64QAM/256QAM| Total 
Q|BFSK 10975 0 0 0 0 0 0 0 0/10975 
 |QFSK 0} 11103 0 0 0 0 0 0 0/11103 
2 \8FSK 0 0} 10982 0 0 119 0 0 0/11101 
= BPSK 0 0 1} 10929 0 0 0 0 0/10930 
8 (QPSK 0 0 0 0} 11062 0 0 0 0/11062 
= BPSK 0 0 0 0 0} 11213 0 0 0/11213 
= [16QAM 0 0 0 0 0 0 11317 0 0/11317 
<= |64QAM 0 0 0 0 0 0 4 9434 1656/1 1094 
256QAM 0 0 0 0 0 0 0 3171 8034|11205 
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Ta 


ble 9. 


Confusion matrix for KPCA classifier in AWGN, SNR = 5 dB. 


Classifier Output (Average Performance = 88.7%, Standard Deviation = 1.0% 








Actual Modulation Type 


Actual Modulation Type 



























































BFSK| QFSK| 8FSK| BPSK) QPSK} 8PSK) 16QAM| 64QAM/256QAM| Total 
BFSK 11182 0 0 0 0 0 0 0 0|11182 
QFSK 0} 11328 0 0 0 0 0 0 0/11328 
8FSK 0 1} 10455 0 0 467 0 0 0|10923 
BPSK 0 0 0} 11078 0 0 0 0 0|11078 
QPSK 0 0 0 0} 11075 0 0 0 0|11075 
8PSK 0 0 0 0 0} 11073 0 0 0|11073 
16QAM 0 0 0 0 0 0 10334 780 11/11125 
64QAM 0 0 0 0 0 0 1431 5872) —3889]11192 
256QAM 0 0 0 0 0 0 493 4199} 6332/11024 

4. GDA 


2 
For GDA, the best kernel was Gaussian: k(x, y) =exp(—2=2| Its overall 


performance at both 20 dB and 5 dB (Tables 10 and 11, respectively) is about the same as 


LDA, but with slightly better performance for 256-QAM and slightly different 


performance for the QAM schemes. 


Table 10. 


Confusion Matrix for GDA classifier in AWGN, SNR = 20 dB. 


Classifier Output (Average Performance = 94.8%, Standard Deviation = 0.6% 


























BFSK|_ QFSK| _8FSK| BPSK] QPSK) 8PSK| 16QAM|/ 64QAM/256QAM| Total 
BFSK 11082 0 0 0 0 0 0 0 0|11082 
QFSK 0} 10981 0 0 0 0 0 0 0|10981 
8FSK 0 0} 10994 0 0 118 0 0 0111112 
BPSK 0 0 0} 11204 0 0 0 0 0/11204 
QPSK 0 0 0 0} 11007 0 0 0 0|11007 
8PSK 0 0 0 0 0} 11238 0 0 0/11238 
16QAM 0 0 0 0 0 0 10677 316 0|10993 
64QAM 0 0 0 0 0 0 0 8830) — 2425)11255 
256QAM 0 0 0 0 0 0 0 2393) 8 735/11128 






































4] 





Actual Modulation Type 


Table 11. 


Confusion Matrix for GDA classifier in AWGN, SNR = 5 cB. 


Classifier Output (Average Performance = 87.7%, Standard Deviation = 1.1% 





























BFSK| QFSK| 8FSK| BPSK) QPSK) 8PSK| 16QAM/ 64QAM/256QAM| Total 
BFSK 11038 0 0 0 0 0 0 0 0|11038 
QFSK 0| 10978 0 0 0 0 0 0 0|10978 
8FSK 0 1} 10617 0 0 412 0 0 0|11030 
BPSK 0 0 0} 11057 0 0 0 0 0|11057 
QPSK 0 0 0 0} 11060 0 0 0 0|11060 
8PSK 0 0 0 0 0} 11120 0 0 011120 
16QAM 0 0 0 0 0 0 8511 2834 15/11360 
64QAM 0 0 0 0 0 0 331 6637, 4293)/11261 
256QAM 0 0 0 0 0 0 85 4331 6680|11096 






































D. PERFORMANCE OF THE PCA AND LDA CLASSIFIERS IN AWGN 
AND FADING CONDITIONS 


Tables 12 through 27 are confusion matrices for the LDA classifier. Tables 28 
through 43 are confusion matrices for the PCA classifier. PCA performs slightly better 


than LDA in most conditions. 

In MATLAB, the functions used to simulate fading effects were rayleighchan.m 
and ricianchan.m. The parameters passed to these functions were a sampling interval of 
1x10°; a maximum Doppler shift of 3.5 Hz for slow fading and 5000 Hz for fast fading; 
path delays of 0 and 1x10” for frequency-flat channels and 0 and 2x10~° for frequency- 
selective channels; average path gains of 0 and —10 dB; and a K-factor of 3 for Ricean 


channels. 
1. LDA 


In a slow, frequency-flat Rayleigh channel, LDA misclassifies some of the PSK 
signals for QAM. The main difference between the results at 20 dB (shown in Table 12) 
and 5 dB (shown in Table 13) is once again additional confusion of the QAM schemes 


for each other. 
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Table 12. 


Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Flat 
Rayleigh Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 86.2%, Standard Deviation = 1.1% 






















































































BFSK| QFSK}  8FSK| BPSK) QPSK) 8PSK}| 16QAM| 64QAM/256QAM| Total 
&|BFSK 11118 0 0 0 0 1 0 0 0/11119 
Zz QFSK 0} 11269 1 0 0 3 0 0 12111285 
2 |BFSK 0 10} 10557 0 0 301 0 0 136]11004 
= \BPSK 0 0 0} 10768 219 8 0 0 010995 
8 (QPSK 0 0 0 0} 10460 0 503 21 149)11133 
= 8PSK 0 0 0 0 0} 10948 0 0 3/10951 
a 16QAM 0 0 0 0 0 0 7173 2169 1838)11180 
<|64QAM 0 0 0 0 0 0 0 3635) —_ 7632/11267 

256QAM 0 0 0 0 0 0 0 782} 10284/11066 
Table 13. Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Flat 
Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 82.3%, Standard Deviation = 1.2% 

BFSK| QFSK} 8FSK| BPSK) QPSK) 8PSK| 16QAM| 64QAM/256QAM| Total 
|BFSK 10940 0 0 0 0 10 0 0 1|10951 
a QFSK 0} 11083 1 0 0 36 0 0 711127 
2 \8FSK 0 9g} 10178 0 0 919 0 0 87|11193 
= \BPSK 0 0 0} 10737 255 13 1 0 0/11006 
8 (QPSK 0 0 0 0} 10393 0 474 22 168]11057 
= PSK 0 0 0 0 0} 11116 0 0 91|11207 
2 16QAM 0 0 0 0 0 15 6369 2683) —2077/11144 
<\64QAM 0 0 0 0 0 17 888 3224, —7051/11180 

256QAM 0 0 0 0 0 21 471 2356} —8287/11135 






































The classifier fares better in Ricean fading conditions due to the presence of the 


line-of-sight path. Once again, most of the difference between the results at 20 dB and 5 


dB (shown in Tables 14 and 15, respectively) is due to additional misclassifications of 


QAM signals. 
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Table 14. 


Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Flat 
Ricean Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 91.1%, Standard Deviation = 0.9% 































































































BFSK| QFSK| 8FSK} BPSK) QPSK) 8PSK} 16QAM) 64QAM/256QAM| Total 

&|BFSK 11208 0 0 0 0 0 0 0 0/11208 
7 |QFSK 0} 11100 1 0 0 10 0 0 0}11111 
2 |BFSK 0 3} 11002 0 0 307 0 0 0/11312 
= BPSK 0 0 0} 10989 57 0 0 0 0/11046 
6 (QPSK 0 0 0 0} 10914 0 135 6 25/1 1080 
= BPSK 0 0 0 0 0} 11064 0 0 0/11064 
a 16QAM 0 0 0 0 0 0 9552 861 564|10977 
<\64QAM 0 0 0 0 0 0 1 5170 5862/1 1033 
256QAM 0 0 0 0 0 0 0 1043) 10126)11169 

Table 15. Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Flat 
Ricean Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 85.5%, Standard Deviation = 1.0% 

BFSK| QFSK|  8FSK} BPSK) QPSK) 8PSK} 16QAM)  64QAM/256QAM| Total 

&|BFSK 11119 0 0 0 0 0 0 0 011119 
 (QFSK O|} 11117 0 0 0 24 0 0 0/11141 
2 \8FSK 0 1} 10347 0 0 765 0 0 111114 
£ BPSK 0 0 0} 11032 56 0 0 0 0/11088 
6 (QPSK 0 0 0 0} 11033 0 141 6 29/11209 
= BPSK 0 0 0 0 0} 11159 0 0 0/11159 
2 16QAM 0 0 0 0 0 1 8181 2059 761|11002 
<= \64QAM 0 0 0 0 0 0 1162 4175 5781111118 
256QAM 0 0 0 0 0 0 561 3133 7356|1 1050 






































In Tables 16 and 17, the classifier’s performance in fast, frequency-flat Rayleigh 


conditions is no better than a guess. The increased Doppler shift causes most signals to 


be misclassified as either 8-PSK or 256-QAM, and there is almost no difference between 


the results at 20 dB (shown in Table 16) and the results at 5 dB (shown in Table 17). 
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Table 16. 


Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleigh Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 11.9%, Standard Deviation = 1.0% 































































































BFSK| QFSK} 8FSK) BPSK) QPSK} 8PSK) 16QAM| 64QAM256QAM| Total 
&|BFSK 0 4 71 0 0 9774 0 0 127011119 
Zz QFSK 0 1 15 0 0 9952 0 0 1207|11175 
2 |BFSK 0 0 4 0 0 9879 0 0 118811071 
= \BPSK 0 0 0 289 160 3799 480 50) 6 190/10968 
8 (QPSK 0 0 0 244 168 3593 581 26, _6430)11042 
= BPSK 0 0 0 296 0} 10610 0 0 134]11040 
a 16QAM 0 0 0} 2574 315 7234 168 91 79011172 
<|64QAM 0 0 0} 2808 243 7376 95 100 601/11223 
256QAM 0 0 0} 2941 195 7355 77 87 535/11190 
Table 17. | Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 11.3%, Standard Deviation = 0.98% 
BFSK| QFSK| 8FSK| BPSK} QPSK) 8PSK| 16QAM| 64QAM/256QAM| Total 
3 BFSK 0 0 67 72 0} 10337 0 0 58411060 
Ee QFSK 0 0 21 74 0} 10555 0 0 595|11245 
2 \8FSK 0 0 1 73 0} 10307 0 0 643]11024 
= \BPSK 0 0 0 320 190 4187 668 49 5697/11111 
8 |QPSK 0 0 0 307 197 4028 668 47| _6059)11306 
= |aPSK 0 0 0 536 0} 10108 2 0 373]11019 
a 16QAM 0 0 0} 3014 419 6666 233 72 602|1 1006 
<164QAM 0 0 0} 4009 336 6376 131 68 335]11255 
256QAM 0 0 0} 4153 281 6048 100 55 337|10974 






































In Tables 18 and 19, the classifier appears to perform much better in Ricean 


conditions even though the fading parameters are otherwise the same as in Tables 16 and 


17. This is because the line-of-sight component preserves the spectral content of the FSK 


signals, resulting in relatively high classification rates for them. The rates for the PSK 


and QAM signals are just as bad as in fast, frequency-flat Rayleigh conditions. 
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Table 18. 


Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat Ricean 


Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 52.6%, Standard Deviation = 1.6% 




























































































BFSK| QFSK| 8FSK) BPSK) QPSK) 8PSK) 16QAM)  64QAM256QAM| Total 

&|BFSK 11014 2 0 0 0 22 0 0 27|11065 
Zz QFSK 0} 10881 5 0 0 124 0 0 131]11141 
2 |BFSK 0 3} 7799 0 0 1600 0 0 1592/10994 
= \BPSK 0 0 o}. 1149 9979 0 5 0 011133 
8 (QPSK 0 0 0 0 34 23 1550 1632 7849/1 1088 
= BPSK 0 0 0 0 0} 10894 0 0 211/11105 
= /16QAM 0 0 0 0 12 161 117 8 10881|11179 
<|64QAM 0 0 0 0 14 255 141 6 10696)11112 
256QAM 0 0 0 0 21 258 147 7, 10750]11183 

Table 19. | Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Flat Ricean 
Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 49.3%, Standard Deviation = 1.5% 

BFSK| QFSK| 8FSK) BPSK) QPSK) 8PSK| 16QAM)  64QAM256QAM| Total 

&|BFSK 11096 1 2 0 0 56 0 0 10/11165 
a QFSK 0} 10669 S) 0 0 314 0 0 3411026 
2 |8FSK 0 1} 5995 0 0 4404 0 0 533/10933 
= \BPSK 0 0 0} 1309 9709 0 0 4 1]11023 
8 (QPSK 0 0 0 0 64 36 2208 1524) 7193/11025 
= 8PSK 0 0 0 3 0} 10227 0 0 897|11127 
a 16QAM 0 0 0 2 47 459 414 24| 1026811214 
<|64QAM 0 0 0 6 56 866 586 S) 9713/11236 
256QAM 0 0 0 2 64 1014 633 14 9524/11251 






































In Tables 20 and 21, the results for slow, frequency-selective Rayleigh conditions 


are shown. At 20 dB, many signals are mistaken for QAM. At 5 dB, 8-FSK is also 


frequently mistaken for 8-PSK. 
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Table 20. | Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Selective 
Rayleigh Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 65.0%, Standard Deviation = 1.4% 

BFSK| QFSK} 8FSK| BPSK); QPSK} 8PSK| 16QAM)  64QAM/256QAM| Total 

o BFSK 10643 5 8 0 0 36 0 0 365|11057 
Zz QFSK 333) 8621 37 0 0 5 0 0} 2302)11298 
2 |BFSK 0 654| 6289 0 0 10 0 0}. 4258/1121 
= BPSK 0 0 0} 9612 1528 27 40 0 40/11247 
8 (QPSK 0 0 0 0 6383 209 2280 207| —-2024/11103 
= BPSK 0 0 0 0 0} 10598 0 0 588/11186 
a 16QAM 0 0 0 0 0 96 1606 1799 7426/10927 
<|64QAM 0 0 0 0 0 60 0 304; 10568)10932 
256QAM 0 0 0 0 0 78 0 38} —10923)]11039 

Table 21. | Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Selective 
Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 63.2%, Standard Deviation = 1.3% 

BFSK| QFSK} 8FSK| BPSK; QPSK} 8PSK| 16QAM)  64QAM/256QAM| Total 

&|BFSK 10916 1 0 0 0 314 0 0 129]11360 
a QFSK 304 8397 22 0 0 1443 0 0 1016)11182 
2 |8FSK 0 573] 5120 0 0 3461 0 0 1851|11005 
= \BPSK 0 0 0} 9257 1548 44 40 3 30|10922 
8 (QPSK 0 0 0 0 6490 188 2269 184 1998)11129 
= 8PSK 0 0 0 0 0} 10094 0 0 877/10971 
2 16QAM 0 0 0 0 1 149 1788 1657; — 7535/11130 
<|64QAM 0 0 0 0 0 182 247 866, —9820/11115 
256QAM 0 0 0 0 0 161 148 596} 10281/11186 






































In Tables 22 and 23, it can be seen once again that the classifier performs better in 


Ricean conditions (all other parameters being the same as in Tables 20 and 21). Many 


signals are misclassified as QAM, and at 5 dB 8-FSK is often misclassified as 8-PSK. 
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Table 22. Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Selective 
Ricean Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 69.3%, Standard Deviation = 1.38% 

BFSK| QFSK| 8FSK} BPSK) QPSK) 8PSK} 16QAM) 64QAM/256QAM| Total 

&|BFSK 10889 1 0 0 0 0 0 0 18811078 
a QFSK 222) 9072 32 0 0 0 0 0 165410980 
2 /8FSK 0 496} 6705 0 0 0 0 0} 3925]11126 
= BPSK 0 0 0} 10194 864 12 28 24 43)11165 
S|QPSK 0 0 0 0 7958 93 1935 224 1018)11228 
= |aPSK 0 0 0 0 0} 10903 0 0 277|11180 
2 16QAM 0 0 0 0 0 56 2037 2052 6814/10959 
<164QAM 0 0 0 0 0 50 0 409} 1063711096 
256QAM 0 0 0 0 0 48 0 37} 11103]11188 

Table 23. Confusion Matrix for LDA Classifier in AWGN Plus Slow, Frequency-Selective 
Ricean Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 67.3%, Standard Deviation = 1.2% 

BFSK| QFSK| 8FSK} BPSK) QPSK) 8PSK} 16QAM) 64QAM/256QAM| Total 

&BFSK 10970 0 1 0 0 125 0 0 51111147 
= QFSK 236) 9065 19 0 0 1047 0 1 892|11260 
2 \8FSK 0 438] 5670 0 0 3076 0 2 1924/11110 
= \BPSK 0 0 0} 10259 817 16 28 29 49)11198 
8 IQPSK 0 0 0 0 7884 96 1855 255 959]11049 
= |BPSK 0 0 0 0 0} 10154 0 0 947/11101 
a 16QAM 0 0 0 0 5 65 2307 1890 6761/11028 
<164QAM 0 0 0 0 1 76 360 1141 9550/11128 
256QAM 0 0 0 0 2 71 202 810 9894/10979 






































The classifier once again breaks down in fast, frequency-selective Rayleigh 


conditions due to Doppler shift, which has a more severe effect than frequency 


selectivity. The results at 5 dB (shown in Table 25) do not differ significantly from those 
at 20 dB (shown in Table 24). 
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Table 24. Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Selective 
Rayleigh Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 11.7%, Standard Deviation = 0.9% 

BFSK| QFSK} 8FSK) BPSK) QPSK} 8PSK) 16QAM| 64QAM256QAM | Total 

&|BFSK 3 1 65 0 0} 10263 0 0 759/11091 
Zz QFSK 0 0 22 0 0} 10298 0 0 698/11018 
2 |BFSK 0 0 3 0 0} 10408 0 0 742111153 
= BPSK 0 0 0 108 128 4252 360 211 6081/11140 
8 (QPSK 0 0 0 90 102 4001 191 125, 6664)11173 
= 8PSK 0 0 0 124 0} 10702 0 0 247|11073 
2 16QAM 0 0 0} 1561 358 8042 119 65 1026)11171 
<|64QAM 0 0 1 1804 262 8115 42 39 794/11057 
256QAM 0 0 0 1955 243 8224 36 22 64411124 

Table 25. Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Selective 
Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 11.2%, Standard Deviation = 1.1% 

BFSK|_ QFSK|_ 8FSK| BPSK| QPSK| 8PSK| 16QAM|__64QAM256QAM| Total 

&|BFSK 2 1 70 62 0} 10589 0 0 517/11241 
a QFSK 0 0 16 71 0} 10408 0 0 483|10978 
2 |8FSK 0 0 6 53 0} 10599 0 0 523/11181 
= \BPSK 0 0 0 154 198 4636 511 289) _5309/11097 
8 (QPSK 0 0 0 148 135 4581 361 238) —5618/11081 
= BPSK 0 0 0 262 0} 10314 0 3 46111040 
a 16QAM 0 0 0} 2290 602 7463 167 62 716/11300 
<|64QAM 0 0 0} 2997 477 7033 108 32 430|11077 
256QAM 0 0 0} 3357 471 6719 79 32 347|11005 






































In Tables 26 and 27, the presence of a line-of-sight component once again assists 


in the classification of FSK signals despite the fast, frequency-selective nature of the 


channel. 
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Table 26. 


Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Selective 
Ricean Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 53.6%, Standard Deviation = 1.6% 
































































































































BFSK| QFSK| 8FSK) BPSK) QPSK) 8PSK) 16QAM)  64QAM256QAM| Total 

&|BFSK 11142 0 0 0 0 20 0 0 36|11198 
Zz QFSK 0} 10956 8 0 0 127 0 0 127|11218 
2 |BFSK 0 3) 7711 0 0 1428 0 0 1917|11059 
= \BPSK 0 0 0} 2635 8432 0 0 0 0/11067; 
8 (QPSK 0 0 0 0 0 0 567 1762) 8831/11160 
= BPSK 0 0 0 0 0} 10093 0 0 1063]11156 
a 16QAM 0 0 0 0 6 8 25 14, 10912/10965 
<|64QAM 0 0 0 0 7 22 35 16) 11048/11128 
256QAM 0 0 0 0 2 20 48 7, 10972]11049 

Table 27. | Confusion Matrix for LDA Classifier in AWGN Plus Fast, Frequency-Selective 
Ricean Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 50.8%, Standard Deviation = 1.6% 

BFSK| QFSK| 8FSK| BPSK) QPSK) 8PSK| 16QAM) 64QAM/256QAM| Total 

&|BFSK 11059 1 2 0 0 39 0 0 4111105 
a QFSK 0} 10831 11 0 0 350 0 0 49)11241 
2 |8FSK 0 2| 5885 0 0} 4527 0 0 658/11072 
= \BPSK 0 0 0} 2804 8214 0 0 0 0/11018 
6 (QPSK 0 0 0 0 53 1 1877 1777 7371|11079 
= 8PSK 0 0 0 0 0} 9800 0 0 1428)11228 
a 16QAM 0 0 0 0 36 129 324 82| 10506)11077 
<|64QAM 0 0 0 0 33 347 478 75) 10111|11044 
256QAM 0 0 0 0 31 435 594 74, 10002}11136 

2. PCA 


PCA performs better than LDA in some cases. In Tables 28 and 29 no significant 


difference can be seen between the results for a slow, frequency-flat Rayleigh channel at 


20 dB and 5 dB, respectively. Most errors are due to the QAM schemes being mistaken 


for each other. 
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Table 28. 


Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Flat 
Rayleigh Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 87.0%, Standard Deviation = 1.0%) 




























































































BFSK| QFSK)  8FSK| BPSK) QPSK) 8PSK) 16QAM)  64QAM/256QAM| Total 

&BFSK 11190 1 1 0 0 5 0 0 011197 
. QFSK 3, 11232 2 0 0 35 0 0 011272 
2 |8FSK 0 17| 10728 0 0 377 0 0 0111122 
= BPSK 0 0 1} 10742 393 1 33 3 13/11186 
8 |QPSK 0 0 0 5| 10052 36 550 62 198]10903 
= aPSK 0 0 0 0 0} 10991 0 0 4110995 
2 16QAM 0 0 0 0 1 103 8219 1372 152511 1220 
<164QAM 0 0 0 0 1 144 310 4570 5962/1 0987 
256QAM 0 0 0 0 2 123 160 1605 9228111118 

Table 29. Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Flat 
Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 86.8%, Standard Deviation = 1.1%) 

BFSK| QFSK)  8FSK| BPSK) QPSK) 8PSK) 16QAM| 64QAM/256QAM| Total 

o| _BFSK| 11051 0 2 0 0 0 0 0 0/11053 
S QFSK 0} 10931 0 0 0 40 0 0 010971 
S 8FSK 0 19} 10884 0 0 415 0 0 011318 
=| BPSK 0 0 1} 10847 352 0 26 4 13]/11243 
3| QPSK 0 0 0 8| 10209 34 572 45 269]11137 
= 8PSK 0 0 0 0 0} 10966 0 0 1|10967 
S 16QAM 0 0 0 0 3 108 8157 1412 1520|1 1200 
%6| 64QAM 0 0 0 0 0 146 319 4577| _5929/10971 
% 256QAM 0 0 0 0 1 152 187 1634 9166/1 1140 






































As with LDA, the PCA classifier achieves better results in Ricean conditions. 
However, the PCA classifier does just as well at 5 dB (shown in Table 31) and 20 dB 
(shown in Table 30). 
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Table 30. 


Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Flat 
Ricean Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 91.5%, Standard Deviation = 1.0%) 




























































































BFSK} QFSK) 8FSK} BPSK) QPSK) 8PSK| 16QAM| 64QAM/256QAM| Total 
&BFSK 11197 0 0 0 0 0 0 0 011197 
. QFSK 2| 11249 0 0 0 21 0 0 0)11272 
2 |8FSK 0 4, 10808 0 0 310 0 0 0111122 
= |BPSK 0 0 1; 11089 85 0 ifs 0 4111186 
6 |QPSK 0 0 0 0} 10664 13 150 13 63/10903 
= 8PSK 0 0 0 0 0} 10995 0 0 0/10995 
2 16QAM 0 0 0 0 0 44 9917 779 480}1 1220 
<164QAM 0 0 0 0 0 42 33 6128 478410987 
256QAM 0 0 0 0 0 44 11 1573 9490111118 

Table 31. Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Flat 

Ricean Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 91.7%, Standard Deviation = 0.9%) 

BFSK} QFSK)  8FSK| BPSK) QPSK) 8PSK| 16QAM| 64QAM/256QAM| Total 
&BFSK 11052 0 1 0 0 0 0 0 0/1 1053 
S QFSK 2| 10959 1 0 0 S) 0 0 0/10971 
2 \BFSK 0 11; 11018 0 0 289 0 0 0111318 
£ BPSK 0 0 Oo} 11141 90 0 8 1 3/11243 
S|QPSK 0 0 0 0} 10925 15 146 15 36/1 1137 
= 8PSK 0 0 0 0 0} 10967 0 0 0/1 0967 
a 16QAM 0 0 0 0 0 33 10030 704 433] 1200 
<\64QAM 0 0 0 0 0 46 33 6089 4803/10971 
256QAM 0 0 0 0 1 40 S) 1605 9485)11140 






































In Tables 32 and 33, the classifier’s performance degrades due to Doppler shift in 


the fast, frequency-flat Rayleigh channel, although, not to the same extent as the LDA 


classifier. Once again, the results are about the same at 20 dB and 5 dB. 
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Table 32. 


Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat 
Rayleigh Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 19.0%, Standard Deviation = 1.2%) 




























































































BFSK| QFSK| 8FSK) BPSK) QPSK) 8PSK)  16QAM| 64QAM/256QAM| Total 
8 IBFSK 1 3 71 0 O} 11122 0 0 0/11197 
= QFSK 0 1 19 0 0} 11252 0 0 0/11272 
2 |8BFSK 0 0 2 0 0} 11120 0 0 0111122 
= \BPSK 0 0 0 0 37 5283 171 91 5604/1 1186 
8 |QPSK 0 0 0 0 37 5295 187 85 5299/10903 
= 8PSK 0 0 0 0 3} 10292 25 14 661/10995 
a 16QAM 0 0 0 5 498 2049 1123 387| — 7158)11220 
<164QAM 0 0 0 13 495 2121 1005 379 6974/1 0987 
256QAM 0 0 0 8 467 2053 1031 407 7152)11118 

Table 33. | Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat 

Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 18.9%, Standard Deviation = 1.2%) 

BFSK| QFSK| 8FSK) BPSK) QPSK)  8PSK) 16QAM|) 64QAM/256QAM| Total 
8 IBFSK 0 0 86 0 0} 10967 0 0 0/11053 
a QFSK 0 0 15 0 0} 10956 0 0 0/10971 
2 \8FSK 0 0 5 0 0} 11313 0 0 0111318 
= \BPSK 0 0 0 0 49 5330 169 88 5607/1 1243 
8 |QPSK 0 0 0 0 57 5407 150 82 5441111137 
= 8PSK 0 0 0 0 8} 10259 26 12 662/10967 
a 16QAM 0 0 0 5 457 2135 1128 425 7050/1 1200 
<|64QAM 0 0 0 13 471 2047 1080 388 6972|10971 
256QAM 0 0 0 17 498 2008 1102 412 7103)11140 






































In fast, frequency-flat Ricean conditions, PCA performs worse than LDA, mainly 


due to additional misclassifications of BPSK as QPSK. Tables 34 and 35 show no 


significant difference between SNRs of 20 dB and 5 dB, respectively. 
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Table 34. | Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat Ricean 
Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 43.7%, Standard Deviation = 1.5%) 

BFSK| QFSK} 8FSK| BPSK) QPSK} 8PSK)| 16QAM|  64QAM/256QAM| Total 

&BFSK 11070 3 0 0 0 50 0 0 0/11123 
. QFSK 0} 10794 6 0 0 219 0 0 011019 
2 |8FSK 0 1} 7844 0 0 3231 0 0 0/11076 
= BPSK 0 0 0 0} 11214 0 24 0 2|11240 
8 |QPSK 0 0 0 0 1 131 762 765 9427/1 1086 
= 8PSK 0 0 0 0 0} 10967 0 0 0/10967 
2 16QAM 0 0 0 0 0 5764 6 3 5456/1 1229 
<164QAM 0 0 0 0 0 7730 5 5 3391/11131 
256QAM 0 0 0 0 1 8079 4 4, 3041111129 

Table 35. Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Flat Ricean 
Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 43.8%, Standard Deviation = 1.7%) 

BFSK| QFSK} 8FSK| BPSK) QPSK} 8PSK)| 16QAM)  64QAM/256QAM| Total 

&BFSK 11020 1 1 0 0 50 0 0 0|11072 
S QFSK 0} 10903 8 0 0 248 0 0 011159 
2 \BFSK 0 0} 7918 0 0 3291 0 0 0/11209 
= \BPSK 0 0 0 0} 11154 0 21 1 111177 
8 |QPSK 0 0 0 0 0 123 731 741 9532|11127 
= 8PSK 0 0 0 0 0} 10942 0 0 0/10942 
a 16QAM 0 0 0 0 0 5647 3 3 5467/11120 
<164QAM 0 0 0 0 0 7675 1 3} 3404/1 1083 
256QAM 0 0 0 0 1 8097 8 3 3002/11111 






































PCA also performs slightly worse than LDA for slow, frequency-selective 


channels, whether a Rayleigh or Ricean channel model is used. In Tables 36 and 37, the 


results for Rayleigh fading are shown at 20 dB and 5 dB, respectively. Once again, the 


performance is nearly identical. 
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Table 36. Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Selective 
Rayleigh Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 62.0%, Standard Deviation = 1.6%) 

BFSK| QFSK) 8FSK| BPSK) QPSK] 8PSK| 16QAM)  64QAM/256QAM| Total 

&|BFSK 10897 3 5 0 0 320 0 0 0/11225 
2 QFSK 302 8476 34 0 0 2138 0 0 0/10950 
2 |BFSK 1 568} 6373 0 0 4369 0 0 011311 
= \BPSK 0 0 1 7889 3050 1 121 7 29/11098 
8 (QPSK 0 0 0 0 5452 241 2655 331 2264/10943 
= PSK 0 0 0 0 0} 11071 0 0 011071 
2 16QAM 0 0 0 0 0 1158 1952 1719 6380}11209 
<|64QAM 0 0 0 0 0 1493 13 585 9104/11195 
256QAM 0 0 0 0 0 1589 6 134) — 9269}10998 

Table 37. Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Selective 
Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 61.9%, Standard Deviation = 1.7%) 

BFSK| QFSK| 8FSK| BPSK) QPSK] 8PSK| 16QAM) 64QAM/256QAM | Total 

&IBFSK 10698 3 2 0 0 337 0 0 0/11040 
= QFSK 325 8564 48 0 0 2202 0 0 011139 
2 \8FSK 0 573} 6340 0 0 4240 0 0 011153 
= \BPSK 0 0 3} 8001 3056 0 106 6 14111186 
8 (QPSK 0 0 0 0 5478 265 2702 355) —-2340)11140 
= PSK 0 0 0 0 0} 11084 0 0 0/11084 
2 16QAM 0 0 0 0 0 1159 1855 1749 6401/11164 
<|64QAM 0 0 0 0 1 1468 11 615 8949/1 1044 
256QAM 0 0 0 0 0 1633 7 137] — 9273}11050 






































selective Ricean channel due to the line-of-sight component. 


In Tables 38 and 39, the classifier performs slightly better in a slow, frequency- 
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Table 38. | Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Selective 
Ricean Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 66.3%, Standard Deviation = 1.5%) 

BFSK| QFSK) 8FSK| BPSK) QPSK]  8PSK) 16QAM) 64QAM/256QAM| Total 

&|BFSK 11043 3 3 0 0 176 0 0 0/11225 
2 QFSK 197 9057 38 0 0 1658 0 0 0/10950 
2 |BFSK 0 471| 6927 0 0 3913 0 0 011311 
= \BPSK 0 0 0} 8934 2112 0 46 0 611098 
8 (QPSK 0 0 0 0 6276 112 2758 277 1520)10943 
= PSK 0 0 0 0 0} 11071 0 0 011071 
a 16QAM 0 0 0 0 0 591 2232 2170 6216]11209 
<|64QAM 0 0 0 0 0 757 1 671 9766/11195 
256QAM 0 0 0 0 0 853 0 104) 10041|10998 

Table 39. Confusion Matrix for PCA Classifier in AWGN Plus Slow, Frequency-Selective 
Ricean Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 66.4%, Standard Deviation = 1.4%) 

BFSK| QFSK| 8FSK| BPSK) QPSK] 8PSK| 16QAM) 64QAM/256QAM| Total 

&IBFSK 10826 1 3 0 0 210 0 0 0/11040 
= QFSK 196 9285 29 0 0 1629 0 0 011139 
2 \8FSK 0 480} 6864 0 0 3809 0 0 011153 
s BPSK 0 0 0 9017 2120 0 40 2 711186 
6 (QPSK 0 0 0 0 6446 126 2827 297 1444111140 
= PSK 0 0 0 0 0} 11084 0 0 0/11084 
a 16QAM 0 0 0 0 0 578 2160 2163 6263/11164 
<|64QAM 0 0 0 0 0 831 0 638 9575/11044 
256QAM 0 0 0 0 0 818 0 116) 10116)11050 






































In Tables 40 and 41, it can be seen that PCA performs better than LDA in a fast, 


frequency-selective Rayleigh channel. Once again, there are no significant differences 


between results at 20 dB and at 5 dB. 
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Table 40. Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Selective 
Rayleigh Fading (SNR = 20 dB). 
Classifier Output (Average Performance = 18.1%, Standard Deviation = 1.3%) 

BFSK| QFSK) 8FSK} BPSK} QPSK}  8PSK| 16QAM|  64QAM/256QAM| Total 
&IBFSK 2 0 61 0 0} 11162 0 0 0/11225 
. QFSK 0 0 12 0 0} 10938 0 0 0/10950 
2 |BFSK 0 0 0 0 O} 11314 0 0 0111311 
= BPSK 0 0 0 0 23 6899 80 36 4060]1 1098 
6 (QPSK 0 0 0 0 19 7561 56 25 3282)10943 
= 8PSK 0 0 0 0 0} 10693 12 4 362/11071 
2 16QAM 0 0 0 3 205 3567 585 244 6605|1 1209 
<|64QAM 0 0 0 2 213 3606 528 213 6633)11195 
256QAM 0 0 1 4 213 3376 575 242 6587|10998 

Table 41. Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Selective 

Rayleigh Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 18.2%, Standard Deviation = 1.4%) 

BFSK| QFSK) 8FSK} BPSK) QPSK) 8PSK|  16QAM|  64QAM/256QAM| Total 
8 IBFSK 0 1 57 0 0} 10982 0 0 0/11040 
I |QFSK 0 0 22 0 O} 11117 0 0 0/11139 
2 \8FSK 0 0 1 0 0} 11152 0 0 0111153 
£ BPSK 0 0 0 1 21 7017 81 44 4022/1 1186 
8 (QPSK 0 0 0 0 19 7609 58 41 3413]11140 
= \gPSK 0 0 0 0 1} 10751 11 8 313/1 1084 
2 16QAM 0 0 0 2 203 3817 547 245 6350/1 1164 
<|64QAM 0 0 0 7 233 3491 578 234 6501|11044 
256QAM 0 0 0 3 223 3354 603 253 6614/11050 






































In Tables 42 and 43, the classifier’s performance is seen to improve due to the 


line-of-sight component in the Ricean channel. The results are not quite as good as those 


for LDA. 
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Table 42. 


Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Selective 
Ricean Fading (SNR = 20 dB). 


Classifier Output (Average Performance = 42.4%, Standard Deviation = 1.5%) 




























































































BFSK| QFSK} 8FSK| BPSK} QPSK} 8PSK)| 16QAM) 64QAM/256QAM| Total 
&BFSK 10998 1 1 0 0 48 0 0 0/11048 
. QFSK 0} 10785 11 0 0 249 0 0 0111045 
2 |8FSK 0 2| 7853 0 0 3464 0 0 0/11319 
= |BPSK 0 0 0 0} 11173 0 2 0 0/11175 
6 |QPSK 0 0 0 0 0 22 153 413] 10402)10990 
= 8PSK 0 0 0 0 0} 10866 0 0 0/10866 
2 16QAM 0 0 0 0 0 5386 0 0 577611162 
<|64QAM 0 0 0 0 0 8653 0 0 241811071 
256QAM 0 0 0 0 0 9432 0 0 1892111324 

Table 43. Confusion Matrix for PCA Classifier in AWGN Plus Fast, Frequency-Selective 

Ricean Fading (SNR = 5 dB). 
Classifier Output (Average Performance = 42.7%, Standard Deviation = 1.6%) 

BFSK| QFSK} 8FSK| BPSK) QPSK} 8PSK) 16QAM)  64QAM/256QAM| Total 
&BFSK 11069 0 2 0 0 52 0 0 0/11123 
S QFSK 0} 10683 5 0 0 258 0 0 0/10946 
2 \BFSK 0 0} 7772 0 0 3357 0 0 0/11129 
= \BPSK 0 0 0 O} 11119 0 2 0 0/11121 
S|QPSK 0 0 0 0 0 22 129 394; 1068411229 
= aPSK 0 0 0 0 0} 11196 0 0 0/11196 
a 16QAM 0 0 0 0 0 5267 0 0 57591 1026 
<\64QAM 0 0 0 0 0 8773 1 0 2356/1 1130 
256QAM 0 0 0 0 0 9090 0 0 2010/11100 






































results for PCA do not vary as much with SNR. 


In summary, LDA performs better than PCA in some channel conditions, but the 


In general, as channel conditions 


deteriorate, the two classifiers tend to mistakenly select either 8-PSK or 256-QAM. 


Ricean channels also generally show better results than Rayleigh channels because the 


line-of-sight component preserves the spectral content of the signal, allowing FSK signals 


to be correctly classified. 


work. 


The next chapter summarizes this thesis and presents recommendations for further 
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VI. CONCLUSIONS 


In this thesis, previous work that investigated the use of higher-ordered moments 
and cumulants in the problem of blind modulation classification was built upon. The 
criteria for selecting moments and cumulants to form feature vectors for the signals were 
discussed. Two classical classification algorithms, PCA and LDA, as well their 
nonlinear, kernel-based equivalents, KPCA and GDA, were applied, and the effects of 
AWGN and multipath fading on resulting performances were investigated. 

Whereas previous work assumed that the statistics would be purely real, this work 
demonstrated the effect of rotating the symbol constellation by a phase offset. This result 
was used to select cumulants whose magnitudes are unaffected by phase shifts. 
Additionally, power normalizations were applied to the cumulants to make them more 
robust to fading effects. 

Nonlinear classifiers turned out not to have any advantage in a problem with low- 
dimension feature vectors. The linear schemes turned out to work better (and with fewer 
computations). The steps taken to make the features more resistant to the effects of 
fading were effective for modest Doppler shift and for frequency-flat channels. The 
classifiers’ performance dropped sharply, however, with increasing Doppler shift. 

The linear classifiers proposed here are very effective in channels that only 
undergo corruption by AWGN and in the sorts of fading channels one might expect from 
a person walking and using a mobile device. On the whole, however, fading remains a 
challenge for blind modulation classification. Future work in this area should focus on 
methods to compensate for fading. 

Additionally, one drawback to the work presented here is that it uses continuous 
streams of symbols. A more sophisticated simulation should attempt to account for 
phenomena above the physical layer of communications systems, such as packetized 
data. Furthermore, simulations in the passband should be conducted to see if this 
problem is tractable when dealing with modern techniques such as spread spectrum 


transmission. 
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APPENDIX A: CUMULANTS EXPRESSED AS FUNCTIONS 


OF EQUAL AND LOWER ORDER MOMENTS 


Table 44. Cumulants Expressed In Terms of Moments. 
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A. 


APPENDIX B: BEHAVIOR OF MOMENTS AND CUMULANTS 


AWGN ONLY 


WITH DECREASING SNR 
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Figure 17. 


E,) in AWGN. 
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Figure 18. 
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in AWGN. 
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Figure 19. E49 in AWGN. 
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Figure 20. E,,, in AWGN. 
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Figure 21. E.,, in AWGN. 
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Figure 22. E.69 11 AWGN. 
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Figure 23. E.;, in AWGN. 
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Figure 24. E.4, in AWGN. 
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Figure 25. E.,; in AWGN. 








E80 








—e— BPSK 
—s— QPSK 
—»— 8PSK 
—<— 16QAM 
—*- 64QAM 
—— 256QAM 














100 40 20 15 12 ~~ & 10 8 7 6 5 3 
SNR [dB] 





Figure 26. E 


x,8,0 


in AWGN. 
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Figure 27. E,,, in AWGN. 
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Figure 28. 
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Figure 29. E.;; in AWGN. 
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Figure 30. E.,, in AWGN. 
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Figure 31. C.49 in AWGN. 
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Figure 32. C,,, in AWGN. 
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Figure 33. 


C.,, in AWGN. 
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Figure 34. 


C,¢9 in AWGN. 
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Figure 35. C,.;, in AWGN. 
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Figure 36. C.4. in AWGN. 
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Figure 37. C,,; in AWGN. 
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Figure 38. C 


x,8,0 


in AWGN. 
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Figure 39. C,,, in AWGN. 
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Figure 40. Cc 


x,6,2 


in AWGN. 
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Figure 41. C,;; in AWGN. 
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Figure 42. C.44, in AWGN. 


AWGN PLUS SLOW, FREQUENCY-FLAT RAYLEIGH FADING 


Parameters for the rayleighchan.m function in MATLAB are: 
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e Sampling interval: 1x10° 
e Maximum Doppler shift: 3.5 Hz 
e Path Delays: [0, 1x10] 


e Average Path Gains: [0, -10] 
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Figure 43. E,,9 1 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 44. E,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 45. 
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ao 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 46. 
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in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 47. E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 48. E69 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 49. E.;, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 50. E..4, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 51. E.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 52. E..39 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 53. E,.,, 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 54. E,4. in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 55. E..;; in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 56. E.,, 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 57. C,.49 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 58. C,.,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 59. C,. in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 60. C6 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 61. C,;, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 62. C,4, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 63. C,.,,; in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 64. C..39 11 AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 65. C,,, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 66. C... in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 67. C,;, in AWGN and Slow, Frequency-Flat Rayleigh Fading. 
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Figure 68. C,44 in AWGN and Slow, Frequency-Flat Rayleigh Fading. 


C. AWGN PLUS SLOW, FREQUENCY-FLAT RICEAN FADING 


Parameters for the ricianchan.m function in MATLAB are: 
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e Sampling interval: 1x10° 

e Maximum Doppler shift: 3.5 Hz 
e = K-factor: 3 

e Path Delays: [0, 1x 107] 


e Average Path Gains: [0, -10] 
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Figure 69. E.,9 1 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 70. E,,, 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 71. E.49 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 72. E,,, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 73. E.,, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 74. E69 in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 75. E.,;, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 76. E 


x,4,2 


in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 77. E,,; in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 78. E.3) 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 79. E.,, 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 80. 


62 in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 81. 


3 In AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 82. E44 in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 83. C,.49 in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 84. C,,, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 85. C,,, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 86. C,6 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 87. C,;, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 88. C,4. in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 89. C,,; in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 90. 


3.9 9 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 91. 
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7, 1n AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 92. C.4, in AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 93. C,;, 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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Figure 94. C,,, 11 AWGN and Slow, Frequency-Flat Ricean Fading. 
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D. AWGN PLUS SLOW, FREQUENCY-SELECTIVE RAYLEIGH FADING 


Parameters for the rayleighchan.m function in MATLAB are: 
e Sampling interval: 1x10° 
e Maximum Doppler shift: 3.5 Hz 
e Path Delays: [0, 2x 10°] 


e Average Path Gains: [0, -10] 
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Figure 95. E59 11 AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 96. E.,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 97. E49) in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 98. E,,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 99. E,,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 100. E69 in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 101. E.;, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 102. 


E.4, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 103. E,,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 104. E..39 in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 105. E.,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 106. E..,. 11 AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 107. E,;; in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 108. E.,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 109. 


C,49 in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 110. 
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3, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 111. 


C,,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 112. 
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x, 


9 in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 113. 


5, 1n AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 114. 
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C, 4, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 115. C.,; in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 116. C..) in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 117. C,,, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 118. C.. in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 119. 


C,;, in AWGN and Slow, Frequency-Selective Rayleigh Fading. 
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Figure 120. 


44 in AWGN and Slow, Frequency-Selective Rayleigh Fading. 


E. AWGN PLUS SLOW, FREQUENCY-SELECTIVE RICEAN FADING 


Parameters for the ricianchan.m function in MATLAB are: 














e Sampling interval: 1x10° 

e Maximum Doppler shift: 3.5 Hz 
e = K-factor: 3 

e Path Delays: [0, 2x 10°] 


e Average Path Gains: [0, -10] 
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Figure 121. E,9 in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 122. E,,, 11 AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 123. E49 12 AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 124. E.,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 125. E..., 11 AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 126. E69 in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 127. E..;, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 128. E..4, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 129. E.,,; in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 130. E... 9 in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 131. E..,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 132. E..42, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 133. E,.;; in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 134. E..,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 135. C,4 ) in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 136. C,,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 137. C,,, in AWGN and Slow, Frequency-Selective Ricean Fading. 


124 





C60 





—e— BPSK 
—s— QPSK 
~*~ 8PSK 
—<— 16QAM 
—*— 64QAM 
—e— 256QAM 














100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 











Figure 138. C,.6 11 AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 139. C.;, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 140. C,.,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 141. C,.,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 142. C..) 1n AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 143. C,.,, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 144. C,6, in AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 145. C,.;; 1 AWGN and Slow, Frequency-Selective Ricean Fading. 
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Figure 146. C.,, 11 AWGN and Slow, Frequency-Selective Ricean Fading. 
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F. AWGN PLUS FAST, FREQUENCY-FLAT RAYLEIGH FADING 


Parameters for the rayleighchan.m function in MATLAB are: 
e Sampling interval: 1x10° 
e Maximum Doppler shift: 5000 Hz 
e Path Delays: [0, 1x10] 


e Average Path Gains: [0, -10] 
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Figure 147. 
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29 in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 148. 
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1, 1n AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 149. 


ao 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 150. 
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in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 151. E,,, in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 152. E69 1 AWGN and Fast, Frequency-Flat Rayleigh Fading. 


x 


132 














E51 





—e— BPSK 
—a— QPSK 
—» 8PSK 
—x— 16QAM 
—*— 64QAM 
—e— 256QAM 














100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 
Figure 153. E.;, in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 154. E..4, 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 155. E,,,; in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 156. 
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in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 157. 


E,,, 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 158. 
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6. in AWGN and Fast, Frequency-Flat Rayleigh Fading. 




















E53 





—e— BPSK 
—s— QPSK 
~*~ 8PSK 
—<— 16QAM 
—*- 64QAM 
—— 256QAM 























100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 
Figure 159. E.;, 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 160. E.,, 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 161. C,.49 1 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 162. C,,, in AWGN and Fast, Frequency-Flat Rayleigh Fading. 


137 























100 


C22 


10 8 


SNR [dB] 








—e— BPSK 
—s— QPSK 
~*~ 8PSK 
—<—— 16QAM 
—*x— 64QAM 
—— 256QAM 








Figure 163. 


22 In AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 164. 
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0 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 165. C,;, in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 166. C.4, in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 167. C,,; in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 168. C..¢9 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 169. C,,, in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 170. C..,. in AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 171. C,;; 11 AWGN and Fast, Frequency-Flat Rayleigh Fading. 
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Figure 172. C,44 in AWGN and Fast, Frequency-Flat Rayleigh Fading. 


G. AWGN PLUS FAST, FREQUENCY-FLAT RICEAN FADING 


Parameters for the ricianchan.m function in MATLAB are: 
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e Sampling interval: 1x10° 

e Maximum Doppler shift: 5000 Hz 
e = K-factor: 3 

e Path Delays: [0, 1x 107] 


e Average Path Gains: [0, -10] 
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Figure 173. E,.,) 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 174. E,,, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 175. E49 119 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 176. E,,;, in AWGN and Fast, Frequency-Flat Ricean Fading. 











E22 








—e— BPSK 
—s— QPSK 
—s— 8PSK 
—»— 16QAM 
—*«— 64QAM 
—e— 256QAM 














100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 
Figure 177. E., in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 178. E69 in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 179. E..;, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 180. E..4, 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 181. E,;, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 182. Eg) 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 183. E 


x,7,1 


in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 184. E64, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 185. E,.;, 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 186. E.,, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 187. C.4) in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 188. C,,, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 189. C.., in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 190. C,¢9 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 191. C,;, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 192. C.4, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 193. C,,; in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 194. C,.. ) 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 195. C,,, 11 AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 196. C.6, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 197. C,;, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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Figure 198. C,,, in AWGN and Fast, Frequency-Flat Ricean Fading. 
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H. AWGN PLUS FAST, FREQUENCY-SELECTIVE RAYLEIGH FADING 


Parameters for the rayleighchan.m function in MATLAB are: 
e Sampling interval: 1x10° 
e Maximum Doppler shift: 5000 Hz 
e Path Delays: [0, 2x 10°] 
e Average Path Gains: [0, -10] 
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Figure 199. E59 11 AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 200. E.,, 11 AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 201. E49 in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 202. E,,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 203. E.,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 204. E69 11 AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 205. E,;, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 206. E.4, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 207. E,,; in AWGN and Fast, Frequency-Selective Rayleigh Fading. 





E80 


yp oOo fF aA DN © 














—e— BPSK 
—s— QPSK 
—»— 8PSK 
—<— 16QAM 
—*- 64QAM 
—*— 256QAM 








0 
100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 
Figure 208. E..39 11 AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 209. E.,, 11 AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 210. E..4, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 211. E.;, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 212. E.,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 213. 49 In AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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C,,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 215. C,,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 216. C,69 in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 217. C,;, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 218. C.,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 219. C,.,; in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 220. C..39 11 AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 221. C,.,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 222. C.,. in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 223. C,;, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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Figure 224. C.,, in AWGN and Fast, Frequency-Selective Rayleigh Fading. 
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I. AWGN PLUS FAST, FREQUENCY-SELECTIVE RICEAN FADING 


Parameters for the ricianchan.m function in MATLAB are: 


e Sampling interval: 1x10° 
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e Maximum Doppler shift: 5000 Hz 
e Path Delays: [0, 2x10° ] 
e Average Path Gains: [0, -10] 
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Figure 225. E,.9 in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 226. E. 
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, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 227. 


E..4) in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 228. 
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in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 229. E,,, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 230. E69 1 AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 231. E.;, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 232. E.,, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 233. E..,, in AWGN and Fast, Frequency-Selective Ricean Fading. 


x, 








E80 








—e— BPSK 
—s— QPSK 
—»— 8PSK 
—<— 16QAM 
—*- 64QAM 
—*— 256QAM 

















100 40 20 #15 12 #10 8 7 6 5 3 
SNR [dB] 
Figure 234. E..39 11 AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 235. E,,, 11 AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 236. E.¢, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 237. E.;, 11 AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 238. E44, 1 AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 239. C,49 1 AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 240. C,,, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 241. C,,, in AWGN and Fast, Frequency-Selective Ricean Fading. 








C60 








—e— BPSK 
—=— QPSK 
~*~ 8PSK 
—— 16QAM 
—*— 64QAM 
—e— 256QAM 














100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 





Figure 242. C,.69 in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 243. C,;, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 244. C.4, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 245. C,,; in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 246. C...9 in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 247. C,,, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 248. C..,. in AWGN and Fast, Frequency-Selective Ricean Fading. 
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Figure 249. C,;, in AWGN and Fast, Frequency-Selective Ricean Fading. 








C44 








—e— BPSK 
—=— QPSK 
~~ 8PSK 
—>— 16QAM 
—*— 64QAM 
—*— 256QAM 














100 40 20 15 12 10 8 7 6 5 3 
SNR [dB] 








Figure 250. C.,, in AWGN and Fast, Frequency-Selective Ricean Fading. 
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APPENDIX C: EFFECT OF CONSTELLATION ROTATION ON 
MOMENTS AND CUMULANTS FOR PSK AND QAM SIGNALS 


A. MOMENTS 


1. PSK Signals 
The complex envelope expression for a PSK signal is given as 
wih) Aeso{ i{ a0 }} (C.1) 


where A is the amplitude, N is the number of symbols, d(k) is a uniformly distributed 
random integer between 0 and N-1. Let x'(k) be the same signal after an arbitrary phase 


shift is introduced: 
x'(k) = sero{ 1[ £404) (C.2) 


The general expression for the (p + q) th-order moment of x'(k) is defined as 


coon fae) on ama] 


which may be rewritten as 


2 as 
Ey. q =XPJG(P -ayare| exp{ i222. a4) 





(C.4) 
= exp( jo(p =G) Ee oes 
2. QAM Signals 
The complex envelope expression for a QAM signal is 
x(k) = a(k) + jb(k), (C.5) 
or, in polar form, 
pat OK) 
k) =,/(a(k)) +(b(k)) exp] j tan! —— |, C.6 
x(k) = \|(a(k))’ + (606) e(4 ) (C.6) 
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where a(k) and b(k) are random variables representing the I and Q channels, 


respectively. Let x'(k) be the same signal after an arbitrary phase shift is introduced: 
x(k) = (alk) +(b(K)) el i{ a" et é}} (C.7) 
a 


The ( pt q) th -ordered moment obtained for this signal is given by 


(owoy +(b(k))’ el i me ‘4 
owoy +(b(k)) exp =i{s H)4))) 
| (a(k)) +(bK)) col {ta san) 


(atk) +(b(k)) ex =) en “))) 


=exp(j6(P-)) Ex ng: 





Xx'P.q 





=exp(j¢(p-q))E (C.8) 


B. CUMULANTS 


From the proof above and the expressions in Appendix A, it will be shown that 


the magnitude of most of the cumulants is unaffected by an arbitrary phase shift. In the 


paragraphs to follow, E. 


~pq tefers to the moment E, ,, whose underlying signal has 


been shifted by an arbitrary phase g. Note that there is no need to consider C,,, and 


C,,, as they are mathematically identical to E,,, and E,,,,. 


x,4,0 


’ je AN 
Cy 40 = E40 ~ aE, ) 


| ee 3 (exp ( j2¢) E09 y 
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2. Co 31 


Cyn = Bai ~3E 202 1 
=exp ( j2o) Eg 3 exp( j2¢) | ea Ieee 
= exp( j2¢)(E, a1 —-3E 20%.) 
= exp(j2¢) Ce 


3. Coo 


Xx 


C,,, 18 not magnitude-invariant to phase shifts. 
4. C60 


Cis0 = Eig —1SE vo can +30(Evoo) 
= exp (j6) E55 —15exp( 72) E, 5 exp ( /40) E49 +30(exp( j20) E,o0) 
= exp( 68) (Eo =15E 3.4L yay +30( E29 )) 
=exp(j69)C, 60 


x,5,1 


Cy51 = E53) LOE,» E31 SE. Ea + 20( E25) Eras 

= exp ( j4¢) | ae 10exp ( j2¢) E99 XP ( j2¢) | a ay Ore exp ( j4¢) Ee 46 
+30 (exp ( j2¢) EL oo y Fea 

= exp ( j4¢) les OE, 4 E31 SE E40 + 30(E,. y Bua 


= €Xp ( j4¢) Cy sae 


C,4» 18 not magnitude-invariant to phase shifts. 


x. 
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Cx 


7. C44 
C..,3; 18 not magnitude-invariant to phase shifts. 
8. C5 


0 = Ego ~35(Exgo) —630(Eino) +420(Ela0) Enso 

= exp( j84) E,,) -35(exp( j4¢) E,,.9) —630(exp(j26)E,o) 
+420(exp(j2¢) Pe) exp(j4¢)E, 4 

= exp( j86)(E,s0 ~35(E,49) —630(E,20). +420(E.20) Exo] 
(780) C, 50. 


9. Coa» Cy629 Cy53> and Cy 44 


C 


er Cx6o» C53, and C,,, are not magnitude-invariant to phase shifts. 
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APPENDIX D: MATLAB CODE 


A. LINEAR CLASSIFIERS 


1 Description: 


e PCA.m performs Principal Component Analysis on training data. 
e LDA.m performs Linear Discriminant Analysis on training data. 


e LINEARTESTER.m generates digitally modulated signals by randomly 
picking the modulation type, then corrupting them with AWGN and 
(optionally) applying fading. Finally, it uses the projection matrix and 
class centroids passed as inputs to classify the signals. Works with either 


class data calculated by either PCA.m or LDA.m 


2. PCA.m 


function [TrainMean centroids projmatrix] = pca(TrainData, Neigvecs, Nclasses, 
Nsamp) 


% Author: Nathan Geisinger, Naval Postgraduate School 
% Date: 24 February 2010 
% PCA.m performs principal component analysis on the column vectors in 


% TrainData. 


% Inputs: TrainData (Data to be used for training the classifier, stored as 
% column vectors) 

% Neigvecs: number of eigenvectors to be used in the projection matrix 
% Nclasses: number of classes present in the training data 


% Nsamp: number of samples per class 


% Calculate mean training vector and subtract from each vector to form a 


% centered training matrix 
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TrainMean = mean(TrainData,2); 


CentTrainData = (TrainData - TrainMean*ones(1,Nclasses*Nsamp)); 


% Calculate covariance matrix and decompose into eigenvectors. Keep the 
% desired number of eigenvectors (Neigvecs) as the projection matrix 
covarmatrix = CentTrainData*conj(CentTrainData)’; 

[eigvecs,eigvals] = eig(covarmatrix); 

SortedEigvals = sort(eigvals,1,'descend’); 

SortedEigvals = SortedEigvals(1,:); 

[SortedEigvals,index] = sort(SortedEigvals,2,'descend’); 


projmatrix = eigvecs(:,index(1:Neigvecs)); 


% Calculate class centroids by projecting training data for each class and 
% calculating the mean 
centroids = zeros(size(projmatrix,2),Nclasses); 
for n = 1:Nclasses 
projclass = projymatrix'*CentTrainData(:,[(n-1)*Nsamp+1:n*Nsamp]); 
centroids(:,n) = mean(projclass,2); 


end 


3. LDA.m 


function [TrainMean centroids W] = Ida(TrainData, Neigvecs, Nclasses, Nsamp) 
% Author: Nathan Geisinger, Naval Postgraduate School 

% Date: 24 February 2010 

% LDA.m performs linear discrimant analysis on the column vectors in 


% TrainData. 


% Inputs: TrainData (Data to be used for training the classifier, stored as 
% column vectors) 
% Neigvecs: number of eigenvectors to be used in the projection matrix 


% Nclasses: number of classes present in the training data 
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% Nsamp: number of samples per class 


% First determine the mean training vector and subtract it from the 

% training data to form the centered training matrix. Initialize the 

% within- and between-class scatter matrices Sw and Sb and the matrix of 
% class centroids. 

r = size(TrainData, 1); 

TrainMean = mean(TrainData,2); 

CentTrainData = (TrainData - TrainMean*ones(1,Nclasses*Nsamp)); 

Sw = zeros(r,r); 

Sb = zeros(r,r); 


centroids = zeros(Neigvecs,Nclasses); 


% Calculate Sw and Sb 
for n = 1:Nclasses 
Si = zeros(r,r); 
mu = mean(CentTrainData(:,(n-1)*Nsamp+1:n*Nsamp),2); 
for m= 1:Nsamp 
x = CentTrainData(:,(n-1)*Nsamp+m); 
Si = Si + (x-mu)*conj((x-mu)’); 
end 
Sb = Sb + Nsamp*(TrainMean - mu)*conj((TrainMean - mu)'); 
Sw = Sw + Si; 


end 


% Decompose inv(Sw)*Sb into eigenvectors and keep the top Neigvecs of them 
% as the projectoin matrix 
[eigvecs, lambda] = eig(inv(Sw)*Sb); 
lambda = sort(lambda, |,'descend'); 
lambda = lambda(1,:); 
[lambda,index] = sort(lambda,2,'descend'); 
W = eigvecs(:,index(1:Neigvecs)); 
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% Find class centroids by projecting each class and taking the mean 
for n = 1:Nclasses 
centroids(:,n) = mean((W'*CentTrainData(:,(n-1)*Nsamp+1:n*Nsamp)),2); 


end 


4. LINEARTESTER.m 


function[tracker,percent]=lineartester(numtrials,msglength,ModSchemes,symorde 


r,SymSNR,MeanProf,centroids,projmatrix,channeltype) 


picking 


% Author: Nathan Geisinger, Naval Postgraduate School 
% Date: 24 February 2010 
% LINEARTESTER.m generates digitally modulated signals by randomly 


% the modulation type, then corrupts it with AWGN and (optionally) fading. 
% Finally, uses the projection matrix and class centroids passed as inputs 


% to classify the signals. Works with either lda.m or pca.m. 


% Inputs: 

% numtrials: number of signals to generate 

% msglength: length of each signal 

% ModSchemes: vector containing number of distinct symbols for each scheme. 
% Usually [2 4 8 2 4 8 16 64 256] for BFSK, QFSK, 8FSK, BPSK, QPSK, 8PSK, 
% 16QAM, 64QAM, 256QAM. 

% symorder: 'gray' or 'binary' 

% SymSNR: SNR used to corrupt signals with AWGN 

% MeanProf: the mean profile of the training data (determined by Ida.m or 

% pca.m), to be subtracted from each signal. 

% centroids: Class centroids determined by Ida.m or pca.m 

% projmatrix: Projection matrix determined by Ida.m or pca.m 


% channeltype: 1 = AWGN only, 2 = Rayleigh, 3 = Ricean, 4 = mixed. 
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NumSchemes = length(ModSchemes); 


% tracker will keep track of the results as a confusion matrix 


tracker = zeros(NumSchemes); 


% define modulator objects for PSK and QAM signals 


bpskmodulator 
modem.pskmod('M',2,'PhaseOffset’,0,'SymbolOrder’,symorder,'InputType'’, integer’); 


qpskmodulator 
modem.pskmod('M',4,'PhaseOffset',0,'SymbolOrder’,symorder, 'InputType'’, integer’); 


epskmodulator 
modem.pskmod('M’,8,'PhaseOffset',0,'SymbolOrder’,symorder, 'InputType'’, integer’); 


modulator 16 
modem.qammod('M',16,'PhaseOffset',0,’SymbolOrder',symorder, InputType'’, integer’); 


modulator64 
modem.qammod('M',64,'PhaseOffset',0,’SymbolOrder',symorder, InputType'’, integer’); 


modulator256 
modem.qammod('M',256,'PhaseOffset',0,'SymbolOrder',symorder,'InputType'’,'integer'); 


% define parameters for FSK signals 


freqsep = 16; 
Nsamp = 16; 
Fs = 128; 


% The following thresholds for distinguishing between BFSK, QFSK, and 8FSK 
% were determined empirically. 

thold23 = 2; 

thold45 = 1.5; 

thold89 = 1.15; 


% Define fading channel parameters 

symrate = 1e6; 

raychan = rayleighchan(1/symrate,5000,[0, 2e-6], [0,-10]); 

richan = ricianchan(1/symrate,5000,3,[0, 2e-6], [0,-10]); 
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h=waitbar(0,'Conducting Trials...'); 
for n = 1:numtrials 
% Pick modulation scheme at random and generate msg 
modtype = randsrc(1,1,1: NumSchemes); 
M = ModSchemes(modtype); 
msg = randsrc(msglength, 1,[0:M-1]); 
switch modtype 
case {1,2,3} 
modmsg = fskmod(msg,M,freqsep,Nsamp,Fs,'cont',symorder); 
case 4 
modmsg = modulate(bpskmodulator,msg); 
case 5 
modmsg = modulate(qpskmodulator,msg); 
case 6 
modmsg = modulate(epskmodulator,msg); 
case 7 
modmsg = modulate(modulator16,msg)./sqrt(10); 
case 8 
modmsg = modulate(modulator64,msg)./sqrt(42); 
case 9 
modmsg = modulate(modulator256,msg)./sqrt(170); 


end 


% Use channel type to determine the propagation model to use 
if channeltype == 4 

propmodel = randsrc(1,1,1:3); 
else 

propmodel = channeltype; 


end 


% Apply the appropriate model. AWGNFreePower is a measurement of the 
% signal power before adding AWGN. 
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switch propmodel 

case 1 
AWGNMFreePower = sqrt(mean(modmsg.*conj(modmsg))); 
noisymsg = awgn(modmsg,SymSNR,'measured’); 

case 2 
multipath = filter(raychan,modmsg(1:msglength)); 
AWGNFreePower = sqrt(mean(multipath.*conj(multipath))); 
noisymsg = awgn(multipath, SymSNR,'measured’); 

case 3 
multipath = filter(@richan,modmsg(1:msglength)); 
AWGNFreePower = sqrt(mean(multipath.*conj(multipath))); 
noisymsg = awgn(multipath, SymSNR,'measured’); 


end 


% Normalize by AWGNFreePower 
noisymsg = noisymsg./AWGNFreePower; 


%First check to see if the signal was FSK 
sfft = fft(noisymsg(1:4096)); 

sortfft = sort(abs(sfft), 1,'descend'); 
ratio23 = sortfft(2)/sortfft(3); 

ratio45 = sortfft(4)/sortfft(5); 

ratio89 = sortfft(8)/sortfft(9); 


if ratio23>thold23 
modguess = 1; 
else if ratio45>thold45 
modguess = 2; 
else if ratio89>thold89 
modguess = 3; 
else 
% Tf not FSK, subtract the mean training vector, project 
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% the result, and calculate the distance to all centroids. 
% Closest distance wins 
sigprof = (profgen(noisymsg) - MeanProf); 
projdata = projmatrix'*sigprof; 
euclidnorm = sqrt(sum((centroids - projdata*ones(1,NumSchemes- 
3)).42)); 
[dist,smodguess] = min(euclidnorm); 
modguess = modguess+3; 
end 
end 
end 
tracker(modtype,modguess) = tracker(modtype,modguess) + 1; 
waitbar(n/numtrials,h); 
end 
close(h); 
% Calculate overall performance 


percent = sum(diag(tracker))/numtrials; 


B. NONLINEAR CLASSIFIERS 


1. Description: 
e KPCA.m performs kernel-based (nonlinear) principal component analysis 
on training data. 
e GDA.m performs general discriminant analysis on training data. 


e NONLINTESTER.m generates digitally modulated signals by randomly 
picking the modulation type, then corrupting them with AWGN and 
(optionally) applying fading. Finally, it uses the paramters calculated by 
either KPCA.m or GDA.m to classify the signal. 
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2. KPCA.m 


function [TrainMean CentTrainData centroids K M1 u] = kpca(TrainData, 
Nclasses, Nsamp, maptype,a,b,d) 


% Author: Nathan Geisinger, Naval Postgraduate School 
% Date: 24 February 2010 
% KPCA.m performs kernel-based (nonlinear) principal component analysis on 


% the column vectors in TrainData. 


% Inputs: TrainData (Data to be used for training the classifier, stored as 
% column vectors) 

% Neigvecs: number of eigenvectors to be used in the projection matrix 
% Nclasses: number of classes present in the training data 

% Nsamp: number of samples per class 

% maptype: Type of kernel. 'gauss'’, 'sig', or 'poly.' 


% a,b,d: kernel paramters (set to 0 if not used by your kernel) 


% Calculate mean training vector and subtract from each vector to form a 
% centered training matrix 
TrainMean = mean(TrainData,2); 


CentTrainData = TrainData - TrainMean*ones(1,Nclasses*Nsamp); 


% Call KERNEL.m to calculate kernel matrix 

K = kernel(CentTrainData,CentTrainData,maptype,a,b,d); 

% M1 = (1/no. of training vectors) * I, used to center the kernelized 
% training data matrix in the higher-dimensional space 

M1 = (1/(Nclasses* Nsamp))*ones(size(K)); 

K = K-M1*K-K*M14+M1*K*M1; 


% Eigen-decomposition of Kernel matrix. Sort vectors and values in 
% descending order and keep any that are at least 1/100 of the magnitude of 
% the largest eigenvalue 


[u,gamma] = eig(K); 
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[sortgamma, index] = sort(diag(gamma),1,'descend'); 

u = u(:,index); 

sortgamma = sortgamma(sortgamma>max(sortgamma)/100); 
gamma = diag(sortgamma); 


u = u(:, 1:length(sortgamma)); 


% K must be reformed from the most significant eigenvectors and eigenvalues 
% in order for the dimensions to match. 

K = u*gamma*u’; 

% Normalize the eigenvectors. This will form the projection matrix 


u = u/(ones(size(u,1),1)*sum(sqrt(u'*K*u))); 


% Find class centroids by projecting the kernel matrix class by class and 
% calculating the mean. 
ProjData = (K*u)’; 
centroids = zeros(size(u,2),Nclasses); 
for n = 1:Nclasses 
projclass = ProjData(:,(n-1)*Nsamp+1:n*Nsamp); 
centroids(:,n) = mean(projclass,2); 


end 


3. GDA.m 


function [TrainMean CentTrainData centroids K M1 alpha] = gda(TrainData, 
Nclasses, Nsamp, maptype, a,b,d) 


% Author: Nathan Geisinger, Naval Postgraduate School 
% Date: 24 February 2010 
% GDA.m performs general discriminant analysis on the column vectors in 


% TrainData. 


% Inputs: TrainData (Data to be used for training the classifier, stored as 
% column vectors) 
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% Neigvecs: number of eigenvectors to be used in the projection matrix 
% Nclasses: number of classes present in the training data 

% Nsamp: number of samples per class 

% maptype: Type of kernel. ‘gauss’, 'sig’, or ‘poly.’ 


% a,b,d: kernel paramters (set to 0 if not used by your kernel) 


% Calculate mean training vector and subtract from each vector to form a 
% centered training matrix 
TrainMean = mean(TrainData,2); 


CentTrainData = TrainData - TrainMean*ones(1,Nclasses*Nsamp); 


% Call KERNEL.m to calculate kernel matrix 

K = kernel(CentTrainData,CentTrainData,maptype,a,b,d); 

% M1 = (1/no. of training vectors) * I, used to center the kernelized 
% training data matrix in the higher-dimensional space 

M1 = (1/(Nclasses* Nsamp))*ones(size(K)); 

K = K-M1*K-K*M14+M1*K*M1; 


% Eigen-decomposition of Kernel matrix. Sort vectors and values in 

% descending order and keep any that are at least 1/100 of the magnitude of 
% the largest eigenvalue 

[u,gamma] = eig(K); 

[sortgamma, index] = sort(diag(gamma),1,'descend'); 

u = u(:,index); 

sortgamma = sortgamma(sortgamma>max(sortgamma)/100); 

gamma = diag(sortgamma); 

u = u(:, 1:length(sortgamma)); 


K = u*gamma*u'; 


% Form the block diagonal matrix W. There will be one block for each class 
% (modulation type). Each block will be equal to (1/no. of samples per 
% class) * ones(no. of samples per class) 
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WI = (1/Nsamp)*ones(Nsamp); 
W=WI; 


for m = 1:Nclasses-1 
W = blkdiag(W,W1); 


end 


% Eigendecomposition of u'Wu. Vectors are used to calculate the projection 
% matrix alpha. 

[beta lambda] = eig(u'*W*u); 

alpha = u*inv(gamma)*beta; 

% Normalize all vectors in alpha. This is the projection matrix. 


alpha = alpha./(ones(size(alpha, 1), 1)*sum(sqrt(alpha'’* K*alpha))); 


% Project kernelized training data and calculate class centroids 
ProjData = (K*alpha)'; 


centroids = zeros(size(alpha,2),Nclasses); 


for n = 1:Nclasses 
projclass = ProjData(:,(n-1)*Nsamp+1:n*Nsamp); 
centroids(:,n) = mean(projclass,2); 


end 


4. NONLINTESTER.m 


function [tracker,percent] = 
nonlintester(numtrials,msglength,ModSchemes,symorder,SymSNR,MeanProf,CentProfs, 
K,M1,centroids,alpha,maptype,a,b,d,channeltype) 


% Author: Nathan Geisinger, Naval Postgraduate School 
% Date: 24 February 2010 
% NONLINTESTER.m generates digitally modulated signals by randomly 
picking 
% the modulation type, then corrupts it with AWGN and (optionally) fading. 
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% Finally, uses the paramters calculated by either KPCA.m or GDA.m to 
% classify the signal. 


% Inputs: 

% numtrials: number of signals to generate 

% msglength: length of each signal 

% ModSchemes: vector containing number of distinct symbols for each scheme. 
% Usually [2 4 8 24 8 16 64 256] for BFSK, QFSK, 8FSK, BPSK, QPSK, 8PSK, 
% 16QAM, 64QAM, 256QAM. 

% symorder: 'gray’ or 'binary' 

% SymSNR: SNR used to corrupt signals with AWGN 

% MeanProf: the mean profile of the training data (determined by Ida.m or 

% pca.m), to be subtracted from each signal. 

% CentProfs: Matrix of centered training profiles (from KPCA.m or GDA.m) 

% K: Kernel matrix (from GDA.m or KPCA.m) 

% M1: matrix used to center the testing signals in the higher-dimensional 

% space (see GDA.m or KPCA.m) 

% centroids: Class centroids determined by GDA.m or KPCA.m 

% alpha: Projection matrix determined by GDA.m or KPCA.m 

% maptype: Type of kernel. ‘gauss’, 'sig', or ‘poly.’ 

% a,b,d: kernel paramters (set to 0 if not used by your kernel) 


% channeltype: 1 = AWGN only, 2 = Rayleigh, 3 = Ricean, 4 = mixed. 


NumSchemes = length(ModSchemes); 
% tracker will keep track of the results as a confusion matrix 


tracker = zeros(length(ModSchemes)); 


% define modulator objects for PSK and QAM signals 


bpskmodulator = 
modem.pskmod('M',2,'PhaseOffset’,0,'SymbolOrder’,symorder,'InputType'’, integer’); 


qpskmodulator = 
modem.pskmod('M',4,'PhaseOffset’,0,'SymbolOrder’,symorder,'InputType'’, integer’); 
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epskmodulator 
modem.pskmod('M',8,'PhaseOffset’,0,'SymbolOrder’,symorder,'InputType'’, integer’); 


modulator16 
modem.qammod('M',16,'PhaseOffset',0,’SymbolOrder',symorder,'InputType'’, integer’); 


modulator64 
modem.qammod('M',64,'PhaseOffset',0,’SymbolOrder',symorder, InputType'’, integer’); 


modulator256 
modem.qammod('M',256,'PhaseOffset',0,'SymbolOrder',symorder,'InputType'’,'integer'); 


% define parameters for FSK signals 


freqsep = 16; 
Nsamp = 16; 
Fs = 128; 


% The following thresholds for distinguishing between BFSK, QFSK, and 8FSK 
% were determined empirically. 

thold23 = 2; 

thold45 = 1.5; 

thold89 = 1.15; 


% Define fading channel parameters 

symrate = 1e6; 

raychan = rayleighchan(1/symrate,3.5,[0, .le-6], [0,-10]); 
richan = ricianchan(1/symrate,3.5,3,[0, .le-6], [0,-10]); 


h=waitbar(0,'Conducting Trials...'); 
for n = 1:numtrials 
% Pick modulation scheme at random and generate msg 
modtype = randsrc(1,1,1: NumSchemes); 
M = ModSchemes(modtype); 
msg = randsrc(msglength, 1,[0:M-1]); 
switch modtype 
case {1,2,3} 
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modmsg = fskmod(msg,M,freqsep,Nsamp,Fs,'cont',symorder); 
case 4 

modmsg = modulate(bpskmodulator,msg); 
case 5 

modmsg = modulate(qpskmodulator,msg); 
case 6 

modmsg = modulate(epskmodulator,msg); 
case 7 

modmsg = modulate(modulator16,msg)./sqrt(10); 
case 8 

modmsg = modulate(modulator64,msg)./sqrt(42); 
case 9 

modmsg = modulate(modulator256,msg)./sqrt(170); 


end 


% Use channel type to determine the propagation model to use 
if channeltype == 4 

propmodel = randsrc(1,1,1:3); 
else 

propmodel = channeltype; 


end 


% Apply the appropriate model. AWGNFreePower is a measurement of the 
% signal power before adding AWGN. 
switch propmodel 
case | 
AWGNFreePower = sqrt(mean(modmsg.*conj(modmsg))); 
noisymsg = awgn(modmsg,SymSNR,'measured’); 
case 2 
multipath = filter(raychan,modmsg(1:msglength)); 
AWGNFreePower = sqrt(mean(multipath.*conj(multipath))); 
noisymsg = awgn(multipath, SymSNR,'measured’); 
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case 3 
multipath = filter(@richan,modmsg(1:msglength)); 
AWGNFreePower = sqrt(mean(multipath.*conj(multipath))); 
noisymsg = awgn(multipath, SymSNR,'measured’); 


end 


% Normalize by AWGNFreePower 
noisymsg = noisymsg./AWGNFreePower; 


%First check to see if the signal was FSK 
sfft = fft(noisymsg(1:4096)); 

sortfft = sort(abs(sfft), 1,'descend'); 
ratio23 = sortfft(2)/sortfft(3); 

ratio45 = sortfft(4)/sortfft(5); 

ratio89 = sortfft(8)/sortfft(9); 


if ratio23>thold23 
modguess = 1; 
else if ratio45>thold45 
modguess = 2; 
else if ratio89>thold89 
modguess = 3; 
else 
% If not FSK, subtract the mean training vector. Calculate 


% the higher-dimensional version of the vector (ksignal) 


% and center in the higher-dimensional space using M1p and 


% M1. 

sigprof = (profgen(noisymsg) - MeanProf); 

ksignal = kernel(sigprof,CentProfs,maptype,a,b,d); 

M1p = 1/(size(ksignal,2))*ones(size(ksignal)); 

ksignal = ksignal - M1 p*K - ksignal*M1 + M1p*K*M1; 
% Project signal and calculate the distance to each 


202 


% centroid. Smallest distance is chosen as the modulation 
% type. 
projsig = (ksignal*alpha)’; 
euclidnorm = sqrt(sum((centroids -  projsig*ones(1,NumSchemes- 
3)).42)); 
[dist,smodguess] = min(euclidnorm); 
modguess = modguess+3; 
end 
end 
end 
tracker(modtype,modguess) = tracker(modtype,modguess) + 1; 
waitbar(n/numtrials,h); 
end 
close(h); 
% Calculate overall performance 


percent = sum(diag(tracker))/numtrials; 


C. SUPPORTING FUNCTIONS 


i. Description 
e PROFGEN.M calculates higher order statistics of the input signal and 
stores them in a vector. 
e TRAINPROFS.m generates PSK and QAM signals for training classifiers. 


e KERNEL.m calculates a kernel matrix given two collection of 


columnwise input vectors. 


e MDIST.m calculates the distance between pairs of vectors stored 


columnwise in two input matrices. 
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2. PROFGEN.m 


function profile = profgen(s) 

% Author: Nathan Geisinger, Naval Postgraduate School 

% Date: 24 February 2010 

% PROFGEN.M calculates higher order statistics of input signal s and stores 
% them in the vector profile. Note that many are commented out since they 


% are not used. 


% First ensure that mean of s is 0 
[r,c] = size(s); 

Ss =8s - ones(r,1)*mean(s); 

sbar = conj(s); 

% spower = mean(s.*sbar); 

% Calculate higher order moments 
e20 = (mean(s.%2)); 

ell = (mean(s.*sbar)); 

e40 = (mean(s.%4)); 

e31 = (mean((s.43).*sbar)); 

%e22 = (mean((s.%2).*(sbar.*2))); 
%e60 = (mean(s.6)); 

e51 = (mean((s.‘5).*sbar)); 

%e42 = (mean((s.*4).*(sbar.*2))); 
%e33 = (mean((s.43).*(sbar.3))); 
e80 = (mean(s.%8)); 

%e71 = (mean((s.%7).*sbar)); 
%e62 = (mean((s.*6).*(sbar.*2))); 
%e53 = (mean((s.%5).*(sbar.3))); 
%e44 = (mean((s.*4).*(sbar.*4))); 


% Calculate cumulants 
%c20 = e20; 
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%cll=ell; 

c40 = e40 - 3*e20.12; 

%c31 = e31 - 3*e20.*el1; 

Teed 2 = 622 62042 = 2 61.142: 

%c60 = e60 - 15*e20.*e40 + 30*e20.43; 

c51 =e51 - 10*e20.*e31 - 5*e11.*e40 + 30*(e20.%2).*el 1; 


Joc42 = e42 - e€20.*e40 - 8¥%el1.%e31 - 6%e20.%e22 + 6%*e20.43 + 
24*(611"2).*e20; 


%c33 = e33 - 6*e20.*e31 - 9*e11.*e22 + 18*(e20.%2).*e11 + 12*e11.43; 
c80 = e80 - 35*e40.42 - 630*e20.%4 + 420*(e20.%2).*e40; 


%c7T1 = e71 - 35*e40.*e31 - 630%*(e20.43).*el1 + 210*e40.*e20.*e11 + 
210*e20.*e31; 


%orO2 = e62 - 15*e40.*e22 - 20*e31.42 + 30*e40.*(e20.%2) + 60*e40.*(e11.42) + 
240*e31.*e11.*e20 + 90*e22.*e20.%2 - 90*e20.%4 - 540*(e20.42).*e1 1.42; 


%ocr53 = e53 - 5*e40.*e31 - 30*e31.*e22 + 90*e31.*e20.42 + 120*e31.*e11.42 + 
180*e22.*e11.*e20 + 30*e40.*e20.*e1 1 -270*(e20.%3).*e11 - 360*(e11.%3).*e20; 


%oc44 = e44 - e40.42 -18*e22.%2 - 16*e31.42 - 54*e20.44 - 144*e11.44 - 


TAZ2* (E2052) Ol LAD: st —12*%e40*020%2° +. 96%631.%e11,%e20 + 
144%e22.*6e11 42 +:72*%e22.*e20.42 +.96%63.1.%e20.%e1 1; 


%c20 = c20./spower; 
%c11 =cl11./spower; 
c40 = (c40).”.5; 
%c31 = (c31).%.5; 
Ted = (622 ),.5¢ 
%c60 = (c60).%(1/3); 
c51 = (c51).4C1/3); 
%oc42 = (c42).%(1/3); 
%c33 = (c33).%(1/3); 
c80 = (c80).4.25; 
%c7T1 = (c71).4.25; 
%c62 = (c62).%.25; 
%ocS3 = (c53).%.25; 


205 


%oc44 = (c44).%.25; 


profile = abs([c40;c51;c80]); 


3. TRAINPROFS.m 


function profmatrix = 
TrainProfs(SymSNR,msglength,symorder,MsgsPerScheme) 


% Author: Nathan Geisinger, Naval Postgraduate School 
% Date: 24 February 2010 
% TRAINPROFS.m generates PSK and QAM signals for training classifiers 


% Inputs: 

% SymSNR: SNR at which the training signals will be generated 

% msglength: length of signals 

% symorder: 'gray' or 'binary' 

% MsgsPerScheme: Number of training signals to generate per modulation 


% scheme 


numstats = length(profgen(1));%This will determine the size of profmatrix 
ModSchemes = [2 4 8 16 64 256];%Number of symbols in each scheme 
NumSchemes = length(ModSchemes); 


% Modem object for each scheme 


bpskmodulator = 
modem.pskmod('M',2,'PhaseOffset',0,'SymbolOrder’,symorder, 'InputType', integer’); 


qpskmodulator = 
modem.pskmod('M',4,'PhaseOffset',0,'SymbolOrder’,symorder,'InputType'’, integer’); 


epskmodulator = 
modem.pskmod('M’,8,'PhaseOffset',0,'SymbolOrder’,symorder, 'InputType'’, integer’); 


modulator 16 = 
modem.qammod('M',16,'PhaseOffset',0,'SymbolOrder',symorder,'InputType'’, integer’); 
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modulator64 
modem.qammod('M',64,'PhaseOffset',0,'SymbolOrder',symorder, InputType'’, integer’); 


modulator256 
modem.qammod('M',256,'PhaseOffset',0,'SymbolOrder',symorder,'InputType'’, integer’); 


h2 = waitbar(0, 'Generating Training Signals...'); 
% Initialize profmatrix, where statistics for each signal will be stored 
profmatrix = zeros(numstats, NumSchemes*MsgsPerScheme); 
for n = 1:length(ModSchemes) 
M = ModSchemes(n); 
for k = 1:MsgsPerScheme 
% Iterate through each modulation scheme, generate signal, add 
% AWGN, and store its statistics in profmatrix 
msg = randsrc(msglength, 1 ,[0:M-1]); 
switch n 
case | 
modmsg = modulate(bpskmodulator,msg); 
case 2 
modmsg = modulate(qpskmodulator,msg); 
case 3 
modmsg = modulate(epskmodulator,msg); 
case 4 
modmsg = modulate(modulator16,msg)./sqrt(10); 
case 5 
modmsg = modulate(modulator64,msg)./sqrt(42); 
case 6 
modmsg = modulate(modulator256,msg)./sqrt(170); 
end 
noisymsg = awgn(modmsg(1:msglength),SymSNR,'measured’); 
profmatrix(:,(n-1)*MsgsPerScheme+k) = profgen(noisymsg); 
waitbar(((n-1)*MsgsPerScheme+k)/(NumSchemes*MsgsPerScheme),h2); 


end 
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end 


close(h2); 


4. KERNEL.m 


function K = kernel(x1,x2,ktype,a,b,d) 

% Author: Nathan Geisinger, Naval Postgraduate School 

% Date: 24 February 2010 

% KERNEL.m calculates the kernel matrix given two collection of columnwise 


% input vectors x1,x2 and kernel parameters. 


switch ktype 
case 'gauss' 
K = exp(-Mdist(x1,x2)./a); 
case 'sig' 
K = tanh(a*x1'*x2+b); 
case 'poly' 
K = (a*x1'*x2+b).4d; 


end 


5. MDIST.m 


function d = Mdist(x1,x2) 

% Author: Nathan Geisinger, Naval Postgraduate School 

% Date: 24 February 2010 

% MDIST.m calculates the distance between columnwise vectors in x1 and x2. 
% In other words, if each vector describes a point, d(m,n) will be the 

% distance between the mth vector of x1 and the nth vector of x2. This was 

% done without for loops using a trick from dist2.m by Christopher Bishop 

% and Ian T Nabney. 


208 


[n,p] = size(x1); 

[m,k] = size(x2); 

xIsq = x1.42; 

x2sq = x2.42; 

sumx 1sq = ones(k,1)*sum(x1sq); 
sumx2sq = ones(p,1)*sum(x2sq); 


d = sumx1sq' + sumx2sq - 2*x1'*x2 


’ 
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